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ABSTRACT 

This report contains the derivation of the differential 
equations of motion of a 3 D.O.F. body joined to a 3 D.O.F. 
body by an elastic tether. The tether is represented by a 
spring and dashpot in parallel. A computer program which 
integrates the equations of motion is also described. 
Although the derivation of the equations of motion are 
for a general system, the computer program is written for 
defining loads in large boosters recovered by parachutes. 
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NOMENCLATURE 

The following is a list of the variables used in the 
computer program with a brief description of each,, 

The notation is displayed in two forms, 1) as it appears 
in the computer program, and 2) as used throughout the 
discussion of this report,, Some of the variables used 
in the report are defined when they are introduced, and 
are therefore not included in the nomenclature „ 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

A - 

A 

Length of tether in X direction from 
decelerator confluence point to fore- 
body confluence point 

m 

AA (6,4) 


Dummy variable used to express incre- 
mental velocities in the numerical 
integration computation 

N /S or 

rad/sec 

AALP (16) 


An array of eight variables representing 



angle-of-attack of the forebody 

rad 

AALPP (16) 


An array of eight variables representing 



angle-of-attack of the decelerator 

rad 

AAM ( 8 ) 


An array of eight variables representing 



Mach number of the forebody 


AAMP ( 8 ) 


An array of eight variables representing 



Mach number of the decelerator 


ABAR 

A 

Distance along longitudinal axis of the 
forebody from the intersection of the 
body axes to the tether confluence point 

9 



positive toward the nose 

m 

AD 

A 

Time derivative of A 

m /sec 

ALP 

a 

Angle-of-attack of forebody 

rad 

ALP DEG 

a 

Angle-of-attack of forebody 

deg 

ALPP 

a 

0 

Angle-of-attack of decelerator 

rad 

ALPPDE 

r' 

a 

P 

Angle-of-attack of decelerator 

deg 
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FORTRAN STANDARD 


ALPPSL 


ALPSL 


AMPSL 


AMSL 


APBAR 


AREAI 


ATMOS 


AOBAR 


DESCRIPTION 

The ratio of two angle-of-attack differ- 
ences used in interpolation of aero- 
dynamic coefficients of the decelerator. 

The ratio of two angle-of-attack dif- 
ferences used in interpolation of aero- 
dynamic coefficients of the forebody. 

Mach number of forebody 

Mach number of decelerator 

The ratio of two Mach number differences 
used in interpolation of aerodynamics 
coefficients of the decelerator 

The ratio of two Mach number differ- 
ences used in interpolation of aero- 
dynamics coefficients of the forebody 

Projection along longitudinal axis of 
decelerator from a line between the 
intersection of body axes and the 
tether confluence point, positive 
toward the nose 

Alphameric input-AREA SEQUENCE OF 
INFLATION 

Alphameric input defining atmosphere 

Projection along longitudinal axis of 
the forebody from a line between the 
intersection of the body axes and the 
bridle confluence point, positive 
toward the nose 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

ALBAR 


Distance along longitudinal axis of the 
forebody from the intersection of the 
body axes to the No. one bridle attach 
point/ positive toward the nose 

m 

A2BAR 

a 2 

Distance along longitudinal axis of the 
forebody from the intersection of the 
body axes to the No. two bridle attach 
point, positive toward the nose 

m 

B 

B 

Length of tether in Z direction from 
decelerator confluence point to fore- 
body confluence point 

m 

BBAR 

B 

Projection along lateral axis of fore- 
body from a line between the inter- 
section of the body axes and the tether 
confluence point, positive up 

m 

BD 

• 

B 

Time derivative of B 

m/sec 

BET1 

01 

Positive angle defined in Figure 6 

rad 

BETlDE 

3 1 

Positive angle defined in Figure 6 

deg 

BET2 

0i 

Positive angle defined in Figure 6 

rad 

BET2DE 

b 2 

Positive angle defined in Figure 6 

deg 

BPBAR 


Projection along lateral axis of decel- 
erator from a line between the inter- 
section of the body axes and the tether 
confluence point, positive up 

m 

BOBAR 


Projection along lateral axis of the fore- 
body from a line between the inter- 
section of the body axes and the bridle 
confluence point, positive up 

m 
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FORTRAN STANDARD 

DESCRIPTION 

UNITS 

BIBAR 


Distance along lateral axis of the fore- 
body from the intersection of the body 
axes to the No. one bridle attach 
point,, positive up 

m 

B2BAR 

B, 

Distance along lateral axis of the fore- 
body from the intersection of the body 
axes to the No. two bridle attach point, 
positive up 

m . 

C 

CA 

C 

C A 

Damping coefficient of tether N 

Axial coefficient of forebody 

sec/m 

CAAP 


Drag area of decelerator 

2 

m 

CAP 

C AP 

Axial coefficient of decelerator 


CCA ( 8 1 16) 


An array of eight by 16 variables 
representing axial force coefficients 
of the forebody corresponding to 
AAM ( 8 ) and AALP (16) 


CCAP f 8 16) 


An array of eight by 16 variables 
representing axial force coefficients 
of the decelerator corresponding to 
AAMP (8 i and AALPP(16) 


CCM (8, 16) 


An array of eight by 16 variables 
representing moment coefficient of the 
forebody corresponsing to AAM (8) and 
AALP (16) 


CCMP (8,16) 


An array of eight by 16 variables 
representing moment coefficients of the 
decelerator corresponding to AAMP (8) and 

AALPP (16) 
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FORTRAN STANDARD 

DESCRIPTION 

•f 

UNITS 

CCMQ (8,16) 


An array of eight by 16 variables 
representing damping moment coefficients 
of the forebody corresponding to AAM(8) 
and AALP (16) 

_ l 

rad" 

CCMQP (8 ,16) 


An array of eight by 16 variables 
representing damping moment coefficients 
of the decelerator corresponding to 
AAMP ( 8 ) and AALPP(16) 

rad"* 1 

CCN (8,16) 


An array of eight by 16 variables 
representing normal force coefficients 
of the forebody corresponding to 
AAM ( 8 ) and AALP (16) 


CCNP (8. 16) 


An array of eight by 16 variables 
representing normal force coefficients 
of the decelerator corresponding to 
AAMP (8) and AALPP(16) 


CHI 

X 

tan A/B 

rad 

CHID 

X 

dX/’dt 

rad/sec 

CHIDDE 

o 

X 

dX/dt 

deg/sec 

CHIDEG 

X 

tan A/B 

deg 

CM 

c 

m 

Moment coefficient of forebody 


CMP 

C 

mp 

Moment coefficient of decelerator 


CMQ 

C 

mq 

Damping moment coefficient of forebody 

- 1 

rad 

CMQP 

C 

mqp 

Damping moment coefficient of decelerator 

- 1 

rad 

CN 

C N 

Normal force coefficient of forebody 
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FORTRA N STANDARD 

CNP C M 

Np 

CONF 


?ANDARD DESCRIPTION UNITS 

Normal force coefficient of decelerator 
Np 

Dummy number used to test for the initial 
configuration of the system 

Dummy number used to test for completion 
of trajectory 

cos 6 cos 6 

cos 0 cos 0 

p p 

d Aerodynamics reference length for forebody m 

dA dA 

d6 d6 m/rad 

dA dA _ A 

)\~ t m/ rad 


CONST 

CTHE 

CTHEP 

D 

DADTHE 

DADTHP 

DAMP 

DBDTHE 

DBDTHP 


DD < 3 , 3 ) 


m/rad 

m/rad 


Distance between the reference center of 
the forebody and the C.g. of the decelerator 

Coefficients of the second derivatives 
in the equations of motion of the 
forebody 

} 

Distance between the two bridle con-< 
nection points on the forebody r 

Distance between the two suspension 
line connection points on the 

decelerator r 


dA 

dA 

d6 

d6 

dA 

dA 

d6 

d3 

P 

P 

<r. 


C • L, 

dB 

dB 

d6 

d0 

dB 

dB 

d'e 

de“ 
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FORTRAN STANDARD DESCRIPTION 


UNITS 


DLS1 

DLS , 

Change in length of first suspension 
line due to tension in the elastic 
system 

m 

DLS2 

DLS 2 

Change in length of second suspension 
line due to tension in the elastic 
system 

m 

DLl 

DL j 

Change in length of the first bridle 
line of the forebody due to 
tension in the elastic system 

m 

DL2 

DL 2 

Change in length of the second 
bridle line of the forebody due 
to tension in the elastic system 

m 

DP 

% 

Aerodynamic reference length of the 
decelerator {same as Ref. dia. D q ) 

m 

DT 

At 

Integration time increment 

sec 

DTI 


Length of inflation time 

sec 

DTP 


Number which controls the number of 
integrations between data output 


DTPC 


Control variable in printout routine 


DTP1 

DTP1 

Input constant which controls the number 
of integrations between data output when 
DT = DTI 

sec 

DTP 2 

DTP 2 

Input constant which controls the number 
of integrations between data output when 
DT ^ DT2 

sec 

DTVC 


Time increment to close thrust valve of 
reaction control system on forebody 

sec 
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FORTRAN 

STANDARD 

DESCRIPTION 


UNITS 

DTI 

DTI 

First integration time increment 


sec 

DT2 

DT2 

Second integration time increment 


sec 

DYPR 

q 

Dynamic pressure acting on forebody 

N/m 2 

DYPRP 
EE (3) 

EPS 


Dynamic pressure acting on decelerator 

Variables representing the forces or 
torque acting on the forebody, shown 
in the equations of motion 

A small positive constant used to check 
for redundant or inconsistent equations 
in CROUT subroutine 

N/m 2 

N or 
m-N 

EPSP1 

£ . 
P‘ 

Positive angle defined in Figure 

6 

rad 

EPSP1D 


Positive angle defined in Figure 

6 

deg 

EPSP2 

£ p/ 

Positive angle defined in Figure 

6 

rad 

EPSP2D 

e p 2 

Positive angle defined in Figure 

6 

deg 

EPS1 

£ ^ 

Positive angle defined in Figure 

6 

rad 

EPS1DE 

£ , 

A 

Positive angle defined in Figure 

6 

deg 

EPS 2 


Positive angle defined in Figure 

6 

rad 

EPS2DE 


Positive angle defined in Figure 

6 

deg 

ETA1 

n. 

Positive angle defined in Figure 

6 

rad 

ETA1DE 

n i 

Positive angle defined in Figure 

6 

deg 

ETA2 

n 

Positive angle defined in Figure 

6 

rad 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

ETA2DE 


Positive angle defined in Figure 6 

deg 

FF (3 ) 


Variables representing the accelerations 
of the decelerator or 

m/s ec 2 
2 

rad/sec 

Q 

g 

Acceleration of gravity of Z 

2 

N/sec 

GAM 

Y 

Flight path angle of forebody 

rad 

GAMDEG 

Y 

Flight path angle of forebody 

deg 

GAMP 

Y P 

Flight path angle of decelerator 

rad 

GAMPDE 

Yp 

Flight path angle of decelerator 

deg 

GR 


Acceleration of gravity at sea level 

m/sec 2 

HHH 


Altitude below which trajectory is ended 

m 

I 


Dummy variable used in DO loops 


IERSW 


Control number used to check for incon- 
sistent or redundant equations in CROUT 
subroutine 


III 


Control variable used in the iteration 
section of SUBR 


HYP (16) 


An array of sixteen variables representing 
the pitch moment of inertia of the decelerator 

2 

corresponding to TTI (16) kg-m 

IY 

J y 

Pitch moment of inertia of forebody 

, 2 
kg-m 

IYP 

x yp 

Pitch moment of inertia of decelerator 

, 2 
kg-m 

J 


Dummy variable used in DO loops 


JJ 


Dummy variable used to control output 


JJJ 


Dummy variable used to control output 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

K 

K 

Spring constant of elastic system 

N/m 

KBL 

KBL 

Longitudinal spring constant of bridle 

N/m 

KBT 

KBT 

Transverse spring constant of bridle 

N/m 

KH1 

KH1 

Spring constant of one bridle line (LH1) 

N/m 

KH2 

KH2 

Spring constant of second bridle line(LH2) N/m 

KKS (8) 

KKS (8) 

Spring constant array of dimension 8 

N/m 

KPHI 

K<f> 

Spring constant of bridle at a pt\ll-off 
angle <£> 

N/m 

KS 

KS 

Spring constant of both decelerator 
suspension lines 

N/m 

KSPKHl 

KSPKHl 

Spring constant of elastic system when 
bridle line 2 is slack 

N/m 

KSPKH2 

KSPKH2 

Spring constant of elastic system when 
bridle line 1 is slack 

N/m 

LAM 

A 

Angular displacement of forebody's 
confluence point using the intersection 
of the forebody's body axes and the 
longitudinal axis as a reference 

rad 

LAMDEG 

A 

Angular displacement of forebody's 
confluence point using the intersection 
of the forebody's body axes and the long- 
itudinal axis as a reference 

deg 

LAMO 


Positive angle defined in Figure 6 

rad 

LAM ODE 


Positive angle defined in Figure 6 

deg 

LAM OP 

A °P 

Positive angle defined in Figure 6 

rad 

LAM'OPD 

A °P 

Positive angle defined in Figure 6 

deg 
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DATE Fe bruary 5 r 19 73 GOODYEAR AEROSPACE 

REV DATE corporation 

■ tilON II. 0»i0 

REV DATE 


FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

LAMl 


Positive angle defined in Figure 6 

rad 

LAMI1DE 


Positive angle defined in Figure 6 

deg 

LAM 2 

X 

2 

Positive angle defined in Figure 6 

rad 

LAM2DE 

A 

2 

Positive angle defined in Figure 6 

deg 

LHl 

LHl 

Length of first bridle line 

m 

LH2 

LH2 

Length of second bridle line 

m 

LSI 

LSI 

Length of first suspension line 

m 

LS2 

LS 2 

Length of second suspension line 

m 

LT 

Lip 

Length of riser line 

dLip 

m 

iVsec 

LTD 

Lp 

dt“ 

LTR 


Length of riser line if bridle is not 




slack 

m 

LT 0 

L m 

To 

Unstretched length of riser line 

m 

LO 

L o 

Distance from intersection of forebody's 
body axes to bridle confluence point 

m 

LOP 

a 

o 

Distance from c.g, of decelerator to 
confluence point of suspension lines 

m 

Ll 

L. 

Distance from intersection of forebody's 
body axis to the negative bridle attach 
point 

m 

L2 

L 

Distance from intersection of forebody's 
body axes to the positive bridle attach 
point 

m 

M 

m 

Mass of forebody 

kg 
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CODE IDENT NO. 


FORTRAN STANDARD 


DESCRIPTION 


UNITS 


MMA (16) 


Added mass of decelerator at T 

An array of 16 variables representing 
added mass of the decelerator corres- 
ponding to TTI(16) 


MUDDEG 


Mass of decelerator at T 


Pull-off angle of riser from suspension 
lines 


MUDEG 


Pull-off angle of riser from suspension 
lines 


NUDEG 


NUPDEG 


PHIB 


Axial g load on forebody (earth g's) 

Axial g load on decelerator (earth g's) 

Normal g load on forebody (earth g's) 

Normal g load on decelerator (earth g's) 

Positive angle defined in Figure 6 

Positive angle defined in Figure 6 

Positive angle defined in Figure 6 

Positive angle defined in Figure 6 

Pull-off angle of riser from forebody's 
confluence point 

Pull-off angle of riser from forebody's 
confluence point used in iteration 
section of SUBR 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

PH1D 

4 > 

dj> 

dt 

rad/sec 

PHIDDE 

• 

<t> 

d(j> 

dt 

deg/sec 

PH1DEG 

<P 

Pull-off angle of riser from forebody's 
confluence point 

deg 

PHI 1 

*i 

Maximum pull-off angle before bridle 
line 2 goes slack 

rad 

PHI IDE 

4 > , 

Maximum pull-off angle before bridle 
line 2 goes slack 

deg 

PHI2 


Maximum pull-off angle before bridle 
line 1 goes slack 

rad 

PHI2DE 

♦a 

Maximum pull-off angle before bridle 
line 1 goes slack 

deg 

POINT (5) 


An array used to transfer data points 
from the program to a tape 


QTHE 

Q e 

Generalized force on 0 equation 

m-N 

QTHEP 

o 

CD 

n 

Generalized force on 0^ equation 

m-N 

QX 

Q x 

Generalized force on X equation 

N 

QXP 

Q xp 

Generalized force on X^ equation 

N 

QZ 

Q z 

Generalized force on Z equation 

N 

QZP 

Q__ 

zp 

Generalized force on Z^ equation 

N 

R 


Radius of planet 

m 

RHO 

P 

Atmospheric density at Z 

kg/m 3 
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CODE IDENT NO. _ 

25500 

FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

S 

S 

Aerodynamic reference area of forebody 

m 2 

SIG1 

a i 

Positive angle defined in Figure 6 

rad 

SIG1DE 

a , 

Positive angle defined in Figure 6 

deg 

SIG2 

a 2 

Positive angle defined in Figure 6 

rad 

SIG2DE 


Positive angle defined in Figure 6 

deg 

SP 

s 

p 

Reference area of decelerator (S ) 

2 

m 

SPI 


Reference area of decelerator during 
inflation 

2 

m 

SSPI (16) 


An array of sixteen variables repre- 
senting reference area of decelerator 
corresponding to TTI(16) 

2 

m 

STHE 

sin 0 

sin 0 


STHEP 

sin 6 

P 

sin 0 

P 


T 

t 

Flight time 

sec 

TC 


Time at which thrust valve on reaction 




control system is closed 

sec 

TDTC 


Time at which DT and DTP change value 




from DTI ■* DT2 and DTP1 ■* DTP2 

sec 

TENS 


Tension in riser line 

N 

THE 

0 

Pitch angle of forebody 

rad 

THED 

0 

d0 

dt 

rad/sec 

THEDDD 

9 » 

0 

d"e 

dt 2 

deg/sec 2 

THEDDE 

0 

d0 

dt 

deg/sec 
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FORTRAN 

THEDEG 

THEDL 

THEDU 

THED2 

THEP 

THEPD 

THPDDD 

THPDDE 


STANDARD 


DESCRIPTION 


e 2 


Pitch angle of forebody 

Forebody’s pitching rate at which 
reaction control thruster begins to 
turn off 

Forebody's pitching rate at which 
reaction control thruster is turned on 

( de _) 2 


P 


Pitch angle of decelerator 

de 

__E 

dt 

d 2 0 
£ 

dt 2 

de 

__E 

dt 


UNITS 

deg 

rad/sec 

rad/sec 

rad 2 /sec 2 

rad 

rad/sec 

deg/sec 2 

deg/sec 


THPDEG 

TI 


e 

p 


TIME I 


TISL 


TOR 


Pitch angle of decelerator deg 

Time at which decelerator inflation 

begins (=0.0) sec 

Alphameric input - TIME SEQUENCE OF 
INFLATION 

Ratio of two time differences used to 
calculate inflation characteristics of 
decelerator 

Maximum value of torque on the forebody 
produced by the reaction control system 


m-N 
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FORTRAN 

STANDARD 

DESCRIPTION 

UNITS 

TORQ 


Value of torque on forebody produced by 
reaction control system at T 

m-N 

TSL 


Ratio of time differences used in the 
interpolation of gust velocity 


TTENS (8) 


Tension array associated with 
KKS (8) 

N 

TTG ( 8 ) 


An array of eight variables representing 
time used in gust interpolation 

sec 

TTI (16) 


An array of 16 variables representing 
time used in the inflation interpolation 

sec 

TTT 


Time at which trajectory is ended 

sec 

V 

V 

Total inertial velocity of forebody 

m/sec 

VD 


Total inertial acceleration of forebody 

m/sec 2 

VG 

v g 

Gust velocity 

m/sec 

VP 

V P 

Total inertial velocity of decelerator 

m/sec 

VPD 


Total inertial acceleration of decelerator 

m /sec 

VS 


Speed of sound at Z 

m/sec 

WG (8) 


An array of eight variables representing 
gust velocity corresponding to TTG (8) 

m/sec 

X 

X 

Horizontal displacement of forebody along 
the ]j£ inertial coordinate 

m 
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FORTRAN STANDARD 


DESCRIPTION 


ZBAR 


Lateral displacement of the c.g. of the 
forebody from the intersection of the 
body axes, positive up 


XPDD 


Hori 2 ontal displacement of decelerator 
along the inertial coordinate 


XBAR 


Vertical displacement of forebody along 
the % inertial coordinate 

Longitudinal displacement of the c.g, of 
the forebody from the intersection of the 
body axes, positive toward the nose 


ZPDD 


Vertical displacement of decelerator 

along the *7 inertial coordinate 
dZ A 

E 

dt 


Ratio of altitude differences used in the 
interpolation of RHO and VS 
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SECTION I - ' INTRODUCTION 


The objective of this report is to present a computer 
simulation of the dynamics of two bodies (coupled by an 
elastic tether) in a plane. This is a simplification of 
a more general problem (see Ref 5). Both bodies have two 
translation degrees of freedom and one rotational degree 
of freedom each; the tether is considered massless and its 
only function is to apply a constraint to the two bodies 
such that they remain in the vicinity of each other. A 
situation in which this simulation would be of use is in 
a deceleration and stabilization study of a re-entry 
vehicle by a parachute system. The re-entry vehicle 
(hereafter denoted as the forebody) is assumed to have 
arbitrary mass and shape characteristics, but the decel- 
erator is considered to be symmetric. Also included in 
the report is a listing and explanation of the computer 
program used to integrate the equations of motion. 
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SECTION II ~ EQUATIONS OF MOTION 

1 o General 

The system is defined as two rigid bodies joined by an 
elastic tether and free to move in a given plane,. Both 
the forabody and the decelerator have 3 D.O.F. In 
general the forebody may have an off center c.g„, but 
the decelerator is considered to be symmetric and homo- 
geneous o The elastic tether is simulated by a spring 
and dashpot in parallel and is attached to the forebody 
by means of a bridle ; the tether is attached to the 
decelerator at the confluence point of the suspension 
lines or at the apex of a BALLUTE. 

The motion is referenced to a Cartesian coordinate 
system fixed on a flat, non-rotating planet,, Coordinate 
system II is an inertial coordinate system (Figure 1); 
XiZ, and X B1 Z fil are body axes for the forebody and 
decelerator respectively „ XZ and X_, Z D are fixed to 
the forebody and decelerator respectively at one point 
and always remain parallel to the inertial 1 Z axes „ 

In general, axes XZ and X 1 Z x intersect at the same point 
but not at the center of gravity of the forebody. There- 
fore, X, Z are not principal axes in general. However, 

i A 

X fi , Z 0i are principal axes. r t is the vector distance 
from the intersection of the body axes (longitudinal 
and lateral axes of the forebody) to the confluence point 
of the forebody. Because of the harness configuration, 
this confluence point changes location discretely or con- 
tinuously during a simulation. This problem will be 
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2 o Kinetic Energy 

Consider an arbitray body rotating and translating in a 
plane (Figure 2)„ Axes ^ are inertial axes; axes 

are orthogonal axes fixed to the body and intersect 
at point 0. Angular velocity, 6, has only one component 
perpendicular to the plane of motion,, Linear velocity 
V 0 has components V 0 , and V 0 1 along the instantaneous 
directions of the XjZj axes respectively,, m is located 
at the center of mass of the body with position X, z" 
relative to the X^j axes (X and Z are considered 
constant in this problem) . u is the velocity of the 
center of mass with respect to the XjZjaxes and it has 
components u , u „ r is the vector from o to m 0 

X Z 1 


z 



FIGURE 2 


RIGID BODY WITH 3 D„ 0 o F. 
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The velocity u can be expressed as follows: 
u = oj x r = (-8 -f)x(X T + Z k ) 


It follows immediately from equation (1): 


u , = -0 Z 

x i 


u , = +0 X 

z i 


If the point o has instantaneous velocity components 

V„ and V along X and Z (V represents a linear 
oxi o z 1 1 10 ^ 

velocity of the body as a whole), the inertial velocity 
along X^Z^are: 


V=V+u = V - 0Z 
xi oxi xi oxi 


V = V + u — V +6x 
Z 1 0 Z 1 Z 1 0 z 1 


The kinetic energy is : 


T = %. m (V 2 + V 2 ) 

xi z 1 


Expanding equation (4); 


^ m [V . 2 + V„ 2 ]+ h m [0 2 (Z 2 + X 2 )] 

' 1 oxi 0 Z 1 


+ m [-V e Z + V 0 X] 

1 0X1 0 Z 1 


m V 0 2 + % I 0 2 + m(-V Z + V X)0 
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Equation (6) applies to both the forebody and the 
decelerator. However, in the case of the deoelerator 
X and Z are zero. Therefore, for the decelerator; 


T = Jjm V 2 + % I 0 2 

p p °p yp P 


Velocities V and V must now be transformed from 

0X1 0 Z 1 

the body axes coordinate system to the inertial 
coordinate system. Figure 3 shows the relationship. 


V A = X cos 0 + Z sin 8 
oxi 


V = -X sin0 + Z cos 8 
o z i 



FIGURE 3 - TRANSFORMATION ANGLE 


The total kinetic energy of the system is: 


T T = h m (X 2 + Z 2 ) + *5 m p (X p 2 + 2 p 2 )+ ^ I y 9 2 + % I yp e p 2 

+ m 0[-X("z cos 0 + X sin Q)-Z(Z sin 0 - X cos 0)] (10) 
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Potential Enerc 


The potential energy in the system is due to the 
weight of the two bodies and to the elasticity in the 
tether. 


J = mg[Z + X sin 8 + Z cos 0]+ m g Z + h K (L T - L Tq ) 2 (11) 

x lr it xx 


L to is the unstretched length of the tether and L T 
is the stretched length of the tether given by the 
geometry of the system. Referring to Figure 1: 


P - P 
l r 2 1 


P a and P 2 are vectors from the inertial coordinate 
system to the confluence points of the forebody and 
decelerator respectively. For the decelerator: 


X i + Z k + r 

p P 2 

A l + B K 
p Pi p pi 

i cos 0 + k sin 8 

P P 


■i sin 9 + k cos 6 

P P 


P, = i (X + A cos 0 - B sin 0 )+ k(Z + A sin 0 + B cos 0 ) 

2 PP PP P PP PPP 


( 17 ) 
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Derivatives of A and B with respect to the coordinates are: 


6A _ 6B 

JT~ 6X 

P 


= 0 


= A sin 0 + B cos 0 


-A sin 0 - B cos 0 

P P P P 


= -A cos 0 + B sin 0 


-nr— = A cos 0 - B sin 0 

60 p p P P P 


Rayleigh's Dissipation Function (See Ref 1 and 2) 

Frictional forces which are proportional to the 
velocity may be derived in terms of a function defined as 


fT= H .Z, C. q. 
** 1=1 i 


where the summation is over all the degree of freedom. 
For this problem, Raleigh damping is considered only 
in the tether. 


J- HcL 
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Lagrange's Equation 

Lagrange's equation for non-conservative forces, 
holonomic, scleronomic constraint, and Rayleigh's 
dissipation function can be written 


(■*£-) - x ^2_ + — = q. (see Ref 1 & 2) (35) 

dt 6^7 6q7 (S q i v 


In equation (35), the term X expresses the 

generalized force exerted by the tether on the "i''th 

6 °I 

degree of freedom. The term is the damping in the 

spring and Q i is the non-conservative aerodynamic and 
reaction control forces „ 

The Lagrangian (L) is equal to the total kinetic energy 
of the system minus the total potential energy of the 
system. 


L “ T T - V T 


Substituting equation (10) and (11) into equation 

(36), the Lagrangian can be written as a function of 

the generalized coordinates, (X, Z, 0, X Z , 0 ). 

P P P 

L = ^ m (X 2 + Z 2 )+ H m (X 2 + Z 2 ) + ^ I 0 2 + k j q ■ 

p p p y yp p 

+ m 0 [ -X ( Z" cos 0 + X sin e) - Z (Z sin 0 - X cos 0 )] 


- mg[Z + X sin 0 + Z cos 0 ] - m p g Z p - % K (L T - L Tfl ) 2 (37) 
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Now operate on equation (37) with equation (35). 


X equation 


= m X - m e(Z cos 0 + X sin e) 


d / fi L «. • § — — 

prr (-T ) = m X - m 0 (Z cos 0 + X sm 0 ) 

oX 


- m 9 2 (-Z sin 0 + X cos 0) 


ia = 

fix 


+ B ff 3 / L 1 


s3T = 

fix 


Z equation 


= m Z - m 0 (Z sin 0 - X cos 0) 


d .fill. 

dt • ' 

5Z 


= m Z - m 0 (Z sin 0 - X cos 0) 


- m 6 2 ( Z cos 0 + X sin 0) 


- - mg 
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( 46 ) 

( 47 ) 

( 48 ) 

( 49 ) 

( 50 ) 

( 51 ) 

( 52 ) 

( 53 ) 

( 54 ) 
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-m g 
P 


5A ^ _ 6B , /T 
[A 6Z + B 6Z ] /L T 


0 Equation 


10 - m[X(Z cos 0 + X sin 0) + Z (Z sin 0- X cos 0)] 


d_ ,6L> 

dt • ' 

r 60 


= 1^0 - m[X(Z cos 0 + X sin 0) 


+ *Z*(Z sin 0 - X cos 0)] + X0 (-Zsin0 + Xcos 0) 


+■ Z 0 (Z cos 0 + X sin 0)] 


=m0[X(-Z sin 0 + X cosO) + Z (Z cos0 + X sin0)] 


-mg [X cos0 - Z sin 0] 


[A ^ A + B ^ 1 / L 
A 60 + B 60 J /L T 


60 


0 


(62) 
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If equations (68) to (71) are substituted into 
equation (35), the resulting equation is: 


K(L t - l T o ) + X + C l t = 0 


X = 


= -[K(L t - l T o ) + C L t ] 


The variable X can now be substituted into equation 
(35) when writing out the differential equations of 
motion„ L T is found by differentiating equation (25) 

i T = ,»* + (74 

A = X + A cos 6 - B sin 0 - X - A cos0 + B sin 0 

P P P P P 


A = X - A 0 sin 0 - B 0 cos 0 - X 

P P P P P P P 


+ A0sin 0 + B0 cos 0 


B = Z + A sin 0 + B cos 0 - Z - A sin0 - B cos0 (77) 

P P P P P 

B = Z + A 0 cos 0 - B 0 sin 0 - Z - A0 cos 0 

P P P P P P P 


+ B 6 sin 0 


Hence X can be expressed as a function of the 
generalized coordinates and their time derivative. 
The six equations of motion are now expressed as 
follows : 
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1) X Equation: 


mk* - m(Z cos 0 + X sin 0)0 = m(-Z sin 0 + X cos 0)0 2 


+ [K(L t - L to ) + C L t ] [j±~] + Q x (79) 

2) Z Equation : 

mZ* - m Cz sin 0 - X cos 0)6* = m (Z cos0 + X sin 0)0 2 - mg 


+ [K(L t - l t0 ) + C L t ] [ £_ ] + Qj 


3) X Equation 


m X 
P P 


= - [K (L„ - L_ ) + C L_] [ ] + Q, 


4) Z Equation 


m p Z p = ' m p 9 - [K(L t - L To ) + C L t ] [ g- ] + Q Zp (82) 


5) 0 Equation: 


I y 0 - m(Z cos 0 + X sin 0) X - m(Z sin 0 - X cos 0)V = 


-mg [X cos 0 - Z sin 0 ]-[k(l t - l T q ) 


+ C L t ] [A 60 + B 6Q ] / L t + Q 0 


6) 0 Equation 


I 0 

yp p 


= “ [K (L m - 


l t.> + c V' A fe 


- + B Q 

p { «p T 8P 
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6. Non-Conservative Generalized Forces 


a) 


Forebody Aerodynamics; 


The aerodynamics of the forebody are given with 
respect to the body axes as shown in Figure 4. 



FIGURE 4 - AERODYNAMICS OF FOREBODY 

Q = -qS (C cos 0 + C sin 0) (85) 

X IN 

Q z = qS (C N cos 0 - C A sin 0) (86) 


The generalized force Q Q is given later in 
equation (90) . 
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l L 
W 
or 



The aerodynamics of the decelerator are given with 
respect to the body axes (Figure 5). 




(87) 

( 88 ) 
(89) 
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c) Reaction Control System : 

A reaction control system may be used to stabilize 
the forebody's pitching motions. This is accomplished 
by checking the pitching rate of the forebody. If 
the absolute value of the rate is above a given upper 
value, a restoring torque (TORQ) is applied to the 
forebody. This restoring torque is maintained 
until a given lower value of pitching rate is reached. 
The torque is then decreased to zero over a finite 
time increment (DTVC). The generalized force Q q is 
now written as: 


= q S d [C + C s (£1)] + TORQ 
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7 o Solution of Equations of Motion 

Because the center of mass of the forebody is not located 
at the intersection of the longitudinal and vertical axes 
(point 0, Figure 2), the equation of motion of the fore- 
body are coupled in the second derivatives. These equations 
((79), (80), (83) have the following form: 


D n X + D u Z + D n 0 = E j 


D 21 X+D 22 Z+D 23 6 = E 2 


D 3i X + D 32 Z + D 33 9 = E 3 


Before numerically integrating equation (91), they are 
separated using Crout reduction (Refer to Ref 3). The 
final form will be: 


q. = f . (x, z, e, x , z , e , x, z, e, x , z , e , t) 

i ' ' ' p' p' p' ' ' ' p' p' p' ' 


i = 1, 2, 3 


The equations of motion for the decelerator are not 
coupled in the second derivative and can be written in 
the form: 

X = F 1 

.? I 

Z = F > 

\ 

0 p = F 3 J 


The six second order differential equations of motion, 
(92) and (93), can now be numerically integrated using 
4th order Runge-Kutta. (Re Ref. 4). 
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SECTION III 

APPLICATION OF THE EQUATIONS OF MOTION TO THE ANALYSIS 
OF A ROCKET BOOSTER RECOVERED BY A PARACHUTE 
(Ref. Figure 6) 

1. General 

The mathematical model defined up to this point applies 
to a general system. Except for the tether line, the 
entire system is rigid. In actual application the 
structure between either body reference point and the 
appropriate tether end is not rigid. In other words there 
is an elastic structure between the end of the tether and 
the referenced body. The tether is attached to the fore- 
body by an elastic bridle, and to the aft body by the 
elastic suspension lines of a parachute. An effective 
system spring constant must be used to adequately account 
for the effect the suspension lines and bridle have on 
the system spring constant. 

The bridle consists of two lines (LH1, LH2) attached to 
points (1 and 2) located on the forebody. The other ends 

of lines LH1 and LH2 attach to the tether. An important 

fact to remember is that the lines can't carry a compres- 
sive load and one bridle line will go slack if the tether 

tension load is directed in such a direction as to lie 
outside c 1 or . Therefore, when the tether tension load 
is directed along one of the bridle lines or outside 
a l or 0 2 ' t ^ ie °PP os; ‘- te goes slack and the tether and 

the line become one longer tether connecting the aft body 
to the forebody at point 1 or 2 depending upon which bridle 
line is carrying the load. The suspension lines in the 
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parachute (aft body) are resolved into two lines (LSI 
and LS2) . Each line has a spring constant (KS) which is 
used together with the bridle spring constant (K<j>, see 
Equation 133) in calculating an effective system spring 
constant (K) for the dynamic two body system. It should 
be noted here that the computer program from which this 
program has been adapted was written for the Viking program. 
This system had a very short rigid tether, the elastic 
effects of which were included in the parachute suspension 
line spring constant (KS). Therefore, one half the tether 
spring constant should be added in series with one of the 
parachute suspension line spring constants, and the result- 
ant spring constant is the spring constant KS used in this 
computer program. 

The parachute shown in Figure 9 can be made elastic or 
rigid by removing or adding the "C" in the comment column 
of the card CALL SUSPEN in the subroutine, SUBR. It has 
been found that a rigid parachute representation results 
in a much faster running program, with little change in 
tether tension when compared to a system with an elastic 
parachute. Therefore, this program calls for the rigid 
parachute simulation; and if desired it can be made elastic 
as discussed above. 
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2. Bridle Geometry 

At any given time the bridle schematic will be as shown 
in Figure 6. The complete geometry of the system, shown in 
Figure 6 can be defined by inputting six variables. These 
variables are LH1, LH2 , Aj, A 2 , B x , B 2 . LH1 , and LH2 are 
always positive and represent the lengths of the two bridle 
lines of the bridle. Aj and A 2 are positive towards the 
nose of the vehicle. They represent the distance to the 
location of the bridle attach points along the centerline. 

Bj and B 2 represent distances to the bridle attach points 
along the lateral axis of the vehicle. With these six values 
the geometry of Figure 6 is defined through the following 
equations: 


■i = /a, 


2 _ 2 
+ B, 


a 2 + B 2 


DL = / (B 2 - B j ) 2 + (A 2 - Aj) 


- 1 2 2 2 / 

! = COS r (L J + DL - L 2 )/(2*L 1 *DL)] 


( 100 ) 
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S 2 


cos" 1 ! (L 2 2 + DL 2 - L j 2 ) / ( 2 *L 2 *DL) ] 

(101) 


= 

cos ' 1 [ (LH 1 2 + DL 2 - LH 2 2 )/ ( 2 *LH 1 *DL) ] 

(102) 

£ 2 

= 

cos -1 [ (LH 2 2 + DL 2 - LH 1 2 )/ ( 2 *LH 2 *DL) ] 

(103) 


= 

TT - 0 j ~ Ej 

(104) 

n 2 

= 

tt 0 2 - e 2 

(105) 

L o 

= 

% 

[ Lj 2 + LH 1 2 - 2 *LHl*L 1 *cos ( 3 j + e x )] 

(106) 

X 2 

= 

tan" 1 ( B 2 /A 2 ) 

(107) 

*1 

= 

tan" 1 (Bj/Aj) 

(108) 


= 

A 2 + cos” 1 [ (L ^ 2 + L Q 2 - LH 2 2 )/( 2 *L 2 *L 0 )] 

(109) 


= 

tt - 6 , - e x - (Aj - X Q ) 

(110) 

a 2 

- 

^ - 0 2 - e 2 - ( A 0 A 2 ) 

(111) 

V 

= 


(112) 


= 

L 0 * COS V 

(113) 

*0 


L 0 * sin v 

( 114 ) 
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Equations (97) to (114) are used to find the point 
which the tension is acting through. If <f> is greater 
than the confluence point of the forebody is located 

at (A i Bj); if cp is less than a 2 , the confluence point 
is at (A 2 , B 2 ) . Otherwise the confluence point is at 
(A 0 , B 0 ). However, (A fl , B 0 ) is a variable depending on 
LH1 and LH2 which depend on the tension and the angle <f> . 
The following method is used to find LH1 and LH2 when each 
leg of the bridle is under tension. Figure 8 shows a 
schematic of the forebody bridle. 

The tension loads in lines LH1 and LH2 are given by 
Equations 115 and 116. 


T j = TENS [ 


sin$ cos a 2 


coscj) sin 02 


sin(a 1 + a ) sin(a x + a 2 ) 


sin^ cosCi cos<|> sina, 

T 2 = TENS [ , — r—] (111 

sin ( ct j+ a 2 ) sin(a x +a 2 ) 

The change in lengths of bridle lines LH1 and LH2 

from unstrained length is given by Equations 117 and 118. 


DL = Tj/KHl 

DL = T 2 /KH2 

The bridle spring constants, KBT and KBL in the 
directions of forces TT and TL are given by Equations 
119 and 120. 


COS O 2 


2 1 


ypni r r ^ . 1 COS O . t. x f 

[{ sin (0 j + <J 2 ) ' KH1 ' ™2 T 


1 W 


sin o 2 2 

KBL = [{ , } 

sin (a x +a 2 ) ; 


1 .sin a, 2 t 

rTriy +{ } * -JL-] 1 

KH1 sin (a x +a 2 ) KH2 J 
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Now using equations 117 and -118 


LHl = LH1 + DL, 


LH2 = LH2 + DL. 


If the value of LHl and LH2 are now used in equations 
(97) to (114) the confluence point will be translated 
and A 0 and B 0 will be used to give 

new values of A, B, <p, y, etc. Now the process is 
repeated. This is done until - *i_i< - 5 degree. 

If the iteration does not converge for i s^lO , the 
program will write out " ITERATION DOES NOT CONVERGE” 
and will continue on. It has been observed that during 
some computer runs the iteration did not converge but 
continued on to the next step without any noticeable 
effect to the results. Usually, if the iteration does 
not converge, a smaller At is needed. This of course 
cost more time on the computer. 
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» Parachute Suspension Geometry 

A typical parachute has many suspension lines. 

To include the effect of each line separately is no 
small task. Consequently, the suspension system is 
assumed to be two lines as shown in Figure 9. 



FIGURE 9 


PARACHUTE GEOMETRY 
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Since the parachute is symmetric, three quantities will 
define its geometry (DLP , LSI, LS2) . The following 
equations result from this input: 


e = cos" [(DLP 2 + LSI 2 - LS2 2 ) /2*DLP*LS1 ) ] (123) 


= cos" [(DLP 2 + LS2 2 - LSI 2 ) /2*DLP*LS2 ) ] 


(124) 


= [ (^-) 2 + LSI 2 - DLP*LSl*cos e . ] ** 

2 n i 


(125) 


cos" 1 [ ( (■— •) + L 2 - LS2 2 )/(DLP*L )] - J (126) 


A = L cos v 

p op p 


(127) 


B = L sin v 
P op P 


(128) 


Like the bridle confluence point the suspension lines 
confluence point can also translate. Summing forces 
in two orthogonal directions, and assuming the system 
to be in equilibrium, yields 


TENS*cos (y + v )= DLS *KS sin e ,+ DLS *KS sin e 
P 1 P 1 2 I 


(129) 


TENS*sin(y + v )=-DLS.*KS cos e + DLS *KS cos e . 

p 1 pi 2 p 2 


( 130 ) 
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I .. 

3 L. 

r w 

4 a: 


Expressing (127) and 128) in matrix form, inverting and 
solving for DLS x and DLS gives: 



Using equation (130) LSI and LS2 are calculated. 


LSI = 

LSI 

+ DLS j 

(133) 

LS2 = 

LS2 

+ DLS 2 

(134) 


If these values are used in equations (123) to (128) the 

confluence point will be translated as shown in Figure 

9, and new values of A and B will be used in the 

P P 

equation of motion. Unlike the bridle it is assumed 
that neither side of the suspension lines will become 
slack. This motion of the suspension line confluence 
point is also included in the iteration process 
mentioned at the end of the preceding Section III-2. 

It is possible to allow each suspension line to stretch 
independently, thereby providing a better simulation. 

To allow the parachute to change geometry under load, 
remove the "C" from the comment column of the card CALL 
SUSPEN in the subroutine SUBR. This allows entry to 
SUSPEN and provides for stretch in the suspension lines. 
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Spring Constant of Elastic System 

When both bridle lines are in tension, the spring 
constant for the bridle is given by: 


Kef) = 


KBL*KBT 


(Ref Pg 28) 


(KBT*cos 2 <J) + KBL*sin 2 <j> ] 

Then the spring constant for the complete system is 

„ _ 2o*KS*K0 . ^ 

K 2 o *KS + K<f> (Ref P 9 22) 

If one bridle line goes slack the spring constant 

becomes either,' 

K = KSPKH1 = 2*KS*KH1 

K KSPKH1 2*KS+KH1 


K = KSPKH2 = 


2 *KS*KH2 
2*KS+KH2 




UJ 

£E 
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SECTION IV - COMPUTER PROGRAM 


1. Inputs 




The format for all numeric inputs is 8F10.0. There 


a couple 

alphameric inputs which use 20A4. The 


following 

is 

a list of all inputs used to make a 


computer 

run 

in the order read in. 


ATMOS 

- 

Alphameric discription of atmosphere. 




1 card 


TIMEI 

— 

Alphameric statement, TIME SEQUENCE 
OF INFLATION, 1 card 


AREAI 

— 

Alphameric statement, AREA SEQUENCE 
OF INFLATION, 1 card 


TTI 


An array of 16 variables representing 
time inflation sequence, 2 cards 

sec 

SSPI 

- 

An array of 16 variables representing 




reference area of decelerator, S , 
corresponding to TTI, 2 cards 

m 2 

MMA 


An array of 16 variables representing 
added mass associated with the decel- 
erator corresponding to TTI, 2 cards 

kg 

HYP 


An array of 16 variables representing 
pitch moment of inertia of the decel- 
erator corresponding to TTI, 2 cards 

kg-m 2 

TTG 


An array of eight variables repre- 
senting time for gust sequence, 

1 card 

sec 

VVG 


An array of eight variables repre- 
senting gust velocities corresponding 
to TTG, 1 card 

m/ sec 
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AAM - An array of eight variables representing 
Mach number of the forebody, 1 card 
AAMP - An array of eight variables representing 
Mach number of the decelerator, 1 card 
AALPE - An array of eight variables representing 

angle-of -attack of the forebody, 1 card deg 
AALPPE- An array of eight variables representing 
angle-of-attack of the decelerator, 

1 card deg 

CCA - An array of eight by eight variables 
representing axial force coefficients 
of the forebody, 8 cards 
CCN - An array of eight by eight variables 

representing normal force coefficients 
of the forebody, 8 cards 
CCM - An array of eight by eight variables 

representing pitch moment coefficients 
of the forebody, 8 cards 
CCMQ - An array of eight by eight variables 

representing pitch damping coefficient 
of the forebody, 8 cards rad 

CCAP - An array of eight by eight variables 
representing axial force coefficient 
of the decelerator, 8 cards 
CCNP - An array of eight by eight variables 

representing normal force coefficients 
of the decelerator, 8 cards 
CCMP - An array of eight by eight variables 

representing pitch moment coefficients 
of the decelerator, 8 cards 



E-ID- 18( 3-70)(J R-2 18) 
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to 
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CCMQP 

- 

An array of eight by eight variables 
representing pitch damping coefficients 
of the decelerator, 8 cards 

_ l 

rad 

TTENS 

•■A 

An array of 8 elements representing 
force in spring KS 

N 

KKS 

- 

An array of 8 elements representing 
the spring constant, KS 

N/m 

X 

- 

Initial range of forebody 

m 

Z 

- 

Initial altitude of forebody 

m 

THE 

- 

Initial pitch angle of forebody 

deg 

THED 

- 

Initial pitching velocity of forebody 

deg/sec 

V 

- 

Initial velocity of forebody 


GAM 

- 

Initial flight path angle of forebody 

deg 

HHH 

- 

Altitude below which trajectory ends 

m 

THEP 

- 

Initial pitch angle of decelerator 

deg 

GAMP 

- 

Initial flight path angle of decelerator 

deg 

VP 

- 

Initial velocity of decelerator 

m/sec 

THEPD 

— 

Initial pitching velocity of 
decelerator 

deg/sec 

TOR 

— 

Maximum value of torque from the 
reaction control system 

m-N 

THEDU 


Forebody's pitching rate at which 
the reaction control thruster is 
turned on giving a torque of TOR 

deg/sec 

THEDL 

“A 

Forebody's pitching rate at which 
the reaction control thruster 
begins to turn off 

deg/sec 
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DTVC 


Length of time for reaction control 




thruster valve to close 

sec 

APBAR 


Distance from c.g. of decelerator 
to confluence point of the 
decelerator suspension lines 

m 

XBAR 

— 

Lateral c.g. off-set of forebody 
positive up 

m 

ZBAR 

— 

Longitudinal c.g. off-set of 
forebody, positive towards nose 

m 

S 

— 

Aerodynamic reference area of 
forebody 

m 2 

D 

— 

Aerodynamic reference length of 
forebody 

m 

M 

- 

Mass of forebody 

kg 

IY 

- 

Pitch moment of inertia of forebody 

kg-m 2 

LTR 

- 

Length of riser line 

m 

C 

— 

Damping coefficient of elastic 
system 

N sec 
m 

DP 

- 

Aerodynamic reference length of 
decelerator (same as Ref. dia. D q ) 

m 

MP 

- 

Mass of decelerator 

kg 

IYP 

- 

Moment of inertia of decelerator 




fully inflated 

kg-m 2 

DTI 

- 

Inflation time 

sec 

T 

- 

Initial time 

sec 

TI 

- 

Time inflation begins 

sec 

DTI 

- 

First integration time increment 

sec 

DT2 

- 

Second integration time increment 

sec 

DTP1 


Number of integrations between 
printout when DT = DTI 
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DTP 2 

Number of integrations between print- 
out when DT=DT2 


TDTC 

Time at which DT1-* DT2 and DTPl->- DTP2 

sec 

TTT 

Time at which trajectory is ended 

sec 

LH1 

Length of first bridle line 



(see Figure 7) 

m 

LH2 

Length of second bridle line 



(see Figure 7) 

m 

AlBAR 

Distance along longitudinal axis 
of the forebody from the intersection 
of the body axes to the negative 
harness attach point, positive toward 



the nose (see Figure 7) 

m 

A2BAR 

Distance along longitudinal axis of 
the forebody from the intersection of 
the body axes to the positive harness 
attach point, positive toward the 



nose (see Figure 7) 

m 

BlBAR 

Distance along lateral axis of the 
forebody from the intersection of 
the body axes to the negative harness 
attach point, positive up 



(see Figure 7) 

m 

B2BAR 

Distance along lateral axis of the 
forebody from the intersection of the 



body axes to the positive harness 



attach point, positive up (See Figure 

7)m 

LSI 

Length of first suspension line 

m 
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LS2 - Length of second suspension line m_ 

DLP - Distance between the two suspension 

line connection points on the 
decelerator m 


KH1 - Spring constant of bridle line LH1 N/m 

KH2 - Spring constant of bridle line LH2 N/m 

HEADER - Title card for plot 
CONT - Condition card for digital program 
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2 . Output 

At predetermined intervals (see inputs eee, fff), the 
following data is outputed. Each term is defined in the 
nomenclature . 

1st row: T, X, XD, XDD , VD, GAMDEG, AM, TENS, PHIDEG, 




QX, 

CA, 

CAP 








2nd 

row: 

TORQ, Z 

, ZD, 

ZDD , 

VPD, 

GAMPDE , 

AMP, 

DAMP, 

PHIDDE , 



QZ, 

CN, 

CNP 








3rd 

row: 

M, 

THEDEG, THEDDE 

, THEDDD , 

NA, 

ALPDEG, DYPR, LTO , 



MUDEG , 

QTHE, 

CM, 

CMP 






4th 

row: 

MP, 

XP, 

XPD , 

XPDD 

, NN, 

ALPPDE , 

DYPRP 

', LT, 

MUDDEG , 



QXP 

, CMQ, CMQP 







5th 

row: 

IY, 

ZP, 

ZPD , 

ZPDD 

, NAP 

, v, 

DP, 

LTD, 

CHIDEG, QZP , 


K, CAAP 

6th row: IYP, THPDEG, THPDDE, THPDDD , NNP, VP, RHO , DCG, 

CHIDDE , QTHEP, C, SPI 

7th ro : LAMDEG, A, B, PHI IDE, PHI2DE, ABAR, BBAR, NUPDEG, 

APBAR, BP BAR, LOP, VG 

When a simulation reaches TTT or HHH , the computer will 
write out" RUN ENDED BY CONSTRAINTS." It will then attempt 
to read in more data cards to initialize for another run. 
If the first card read in contains a "1." in the first two 
columns, the program will CALL EXIT. Otherwise it will 

I read in data starting from input (y) and proceeding to 

input (www) . 
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Before beginning each trajectory, a list of variables 
will be printed out which mostly define the initial 
geometry of the system. These variables are defined 
in the nomenclature. 

1st row: LHl , LH2 , A1BAR, A2BAR, BIBAR, B2BAR, Ll, L2, 

DL, L0, BETIDE, BET2DE 

2nd row: EPS1DE, EPS2DE, ETA1DE , ETA2DE, SIG1DE, SIG2DE, 

LAMIDE, LAM2DE, LAMODE, NUDEG, AOBAR, BOBAR 

3rd row: K, KSPKH1, KSPKH2 , TI , DTI, THEDU , THEDL , TOR, 

DTVC , LSI, LS2 , DLP 

4th row: LOP, LAMOPD, NUPDEG, APBAR, BP BAR, EPSP1D, 

EPSP2D, DT, KH1 , KH2 

Also printed out is the atmosphere used; the area versus 
time inflation sequence; spring constant array, KKS(8), 
and its tension array, TTEN S ( 8 ) ; and the aerodynamic 
coefficient arrays and their associated mach number and 
angle of attack arrays. 
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Fortran IV Program Description 


A listing of the program may be found in Section IV- 5. 

A description of the main program and the subroutines 
follows : 

a. Main Program 

1) Read inputs 

2) Initialize and define certain variables 

3) Call BRIDLE so that the geometry of the system 
may be defined 

4) Output data points on tape for later plotting if 
T > TPLOT 

5) Check time and altitude constraints (inputs ee 
and hhh) 

6) If constraint or output conditions are met, 
write output 

7) Advance the six coordinates through one time 
increment (At) by use of Runge-Kutta numerical 
integration 

8) Check for increase in DT and DTP (input ggg) 
b o Subroutine SUBR 

1) Calculate acceleration of gravity at Z 

2) Calculate torque available from reaction control 
system 

3) Calculate gust velocity as a function of time, 
and density and speed of sound as a function of 
altitude 

4) Calculate total inertial velocity, Mach number, 
and dynamic pressure of the forebody and 
decelerator 




EF: ENGINEERING PROCEDURE S-017 


PAGE 43 

gee- 15853 

CODE IDENT NO. 


date February 5, 1973 

REV DATE 

REV DATE 


GOODYEAR AEROSPACE 

CORPORATION 

uion IS. OHIO 



5) Calculate flight path angle and angle-of-attack 

6) Call AERO to determine the aerodynamics of the 
system 

7) Calculate bridle, riser, suspension lines geometry 
under tension and iterate until the geometry 
converges to one compatible to the tension in the 
elastic system. After the iteration has converged, 
the tension, damping, pull-off angles and rates 
will be determined. 

8) Call MATRIX to determine the equation of motion of 
the forebody 

9) Express equations of motion of decelerator 

c) Subroutine AERO 

1) Calculate the aerodynamic coefficients of the 
forebody and the decelerator as a function of 
Mach number and angle-of-attack 

2) Calculate the generalized forces acting on the 
forebody 

3) Calculate the reference area, apparent mass, pitch 
moment of inertia and total mass of decelerator 

4) Calculate the aerodynamic coefficients of the 
decelerator (CNP, CMP, CMQP) during inflation, 
assumed a linear increase 

5) Calculate the generalized forces acting on the 
decelerator 

d) Subroutine MATRIX 

1) Define the matrices DD(i,j) and EE(i) 

2) Call CROUT so that the equations of motion of the 
forebody are separated into a form suitable for 
integration 
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e) Subroutine BRIDLE 

1) Calculate the geometry of the bridle, bridle 
attach points and intersection of the body axes 
of the forebody 

2) Call SUSPEN so that the geometry of the decelerator 
and suspension lines can be defined. 

3) Write out pertinent information about the geometry 
of the system 

4) Define the initial position of the decelerator 
with respect to the forebody 

f ) Subroutine SUSPEN 

1) Calculate the geometry of the decelerator and 
suspension lines 

2) Calculate the position of the decelerator if 
the suspension lines are allowed to stretch 

g ) Subroutine CROUT 

1) Decouple the equations of motion of the forebody. 
Equations (91) will be reduced to the form (92) 


4. 


S ampl e C omputer Run 

Figure 10 shows Calcomp plots of 8 program output variables 
versus time for a sample computer run. Immediately 
following this figure is a list of the input data used to 
make this run and several pages of output. 
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0.0 

0 . 0 

0.0 

0.0 

0.0 

0.3 

. 0.0 

0 - 0 . 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0 . 0 

Q . q 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 . 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0.0 

0.0 

0.0 



AFTBODY 

NORMAL 

COEF. 

ARRAY t_ 

CCNP ( 8 i 1 6 ) 










0.0 

0.0690 

0.6200 

0 .0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.3 

0.0 

0.0690 

0.6200 

0.0 

0 . 0 

0 . 0 

0.0 

o.o 

0.0 

0 . 0 

0.0 

0.0 

O.o 

0.0 

0.0 

o.o 

0.0 

0.0 . .. 

0.0 ... 

0 . 0 

0 . o _ 

0.0 

Q.Q 

0.0 

0.0 

0.0 

o_.o 

0,0 

0.0 

0.0 

0.0 

0 . 0 

0 . 0 

0 . 0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0 . 0 

0.0 

0.0 

0.0 

0 .0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0.0 

o.o 

0.0 

0 . 0 

JL-9 

JLiQ 

Q.Q 

0.0 

0.0 

_Q..Q 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

Q.Q 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0.0 

_ 0 . 0 _ 

0.0 
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AFT BO DY P I TCH M OM COEF. ARRAY, CCMP{ 8, 16) 


0 . 0 

- 0.1000 - 

0.9000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0.0 

0.0 

0.0 

0.0 

- 0.1000 - 

0.9000 

0.0 

0.0 

0.0 

O.Q 

0,0 

0.0 

0.0 

0.0 

0.0 

Q , 0 

0.0 

0.0 

JL .. Q _ 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0.0 

0 . 0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0 . 0 

0 . 0 

o . o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0.0 . . 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

... 0.0 ... 

0.0 

0 . 0 

0 . 0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

_ 9 .Q 

0 . 0 
il.q 

0 . 0 
JL .0 

0 . 0 

0.0 

0 . 0 
_Q 

0.0 
_ Q . 0 . 


AFTBODY P I TCH. OAM PI NG COEF , ARRAY, CCMQP(8 , 16) 


- 0 . 1000 

.- 0.10 00 

- 0. 1000 

0.0 

0.0 

0.0 

- _ 0 . 0 __ 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0 . o 

_. 0.0 

0.0 

- o . 1000 

_ ro,iaoo 

- 0.1000 0 .0 

0.0 

0.0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0-0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 * 

0.0 

0.0 

0.0 

0.0 

o . 6 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

... 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

_ 0.0 

_ ...0 . 0 

0.0 

0.0 . 

0.0 

0.0 

o-o 

0.0 

_ 0 . 0 . 

— 0.0 

0 . 0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0-0 .... 

0.0 
0 - 0 ... 
0.0 
0 . 0 

a.Q 

0.0 __ 


0.0 

. 0.0 

0.0 

_ 0 . 0 

_9..Q 







TIME RANGE _ HORVEL HO R ACC TOT ACC GAMMA MACH_.NO TE NS I 0N__._ PHI QX CA CAP 

TORQUE AL T I TUDE.__ VER TVEL VE RT ACC __ . TOT AC.CP GAMMAP MACH NOP__. DAMPING . S HID Q2 CN CN° 

MASS THETA THEVEL THE ACC AX-G ALPHA DY-PR. UTQ MU QT HE CM CMP 

MASSP RANGE P HOKVELP. _HORACCP NOrt-G ALPHAP OY-PRP LT MUD QXP CMQ. CMQP._ _ 

IV ALTI TUDEP. VER TVELP .VERT AC C£ AX— GP VELOCITY. _QP .. LTD CHI ..I. QZP , K ”... CAAP__~ 

2Y_ P I H ETA P T HE V£LP_ THE'aCC P N Ok-GP _V E LQCUY£_DL N SITY ~ D C G CH I 0 01 H FP f. SP 

L AMD A. A R ... PH u PHI2 ABAR RHAR NUP ... A PsJ AR DPDAR LOP GUST _.YEL 


0 . 0 

0,0 

65.3 

.- 6.9 

9 . 9 ... 

- 76.0 . 

0.6 3685.7 

- 6.6 

- 661099 . 

- 0. 72111 

0.50000 

0.0 

6096.0 

- 181 . 9 

... - 8 , 6 

4.4 

- 76.0 

0.6 3.3 

0.0 

... 102263 ... 

3 . 1111 - 

- 0.0000 

9 0682 .20 

90.00 

0 .0 

o 

a 

- 0.38 

166.00 11676.61 15.26 

Q . 0 3 

22 7216 11.6366 

- 0.0033 

2608. 055 

- 28. 22 

65.35 

rl . 26 __ 

0.50 

- 0.03 

11676 . 61 . 15.25 

0.0 

- 3868 . 

- 0 . 0 

- 0 . 1000 

26606760 . 

. 6200.8 

- 13 1.9 

- 6.2 

_ 0.6 

187 . 5 

66.0 0.0 

- 16.0 

15515 . 

557480 . 

1.4 

566732. 76 

- 76. 00 

0.0 

0.00 _ 

- 0.23 

187.65 

0 ,653 ID 00 108.50 

.. 0.0 

- 0 . _ 

875.63 

2.79 

18-66 

- 3.69 

16 . 79 

9 . 08 

76.10 

23.88 . 

8 . 06 0.00 

68.12 

3.00 

. 68.12 

0.0 ... 

,_ Q . 050 _ 


65 . 1 

- 6 . 9 

9 . a 

- 76.1 

0.6 9383.1 

- 6.7 

-639073 . 

- 0.71983 

0.60000 

0.0 

6086.9 

- 182.3 

- 8.5 

15.5 

- 76.0 

0.6 589.1 

. - 3.9 

102568 . 

.. 3.0379 

0.0332 

90682.20 

90.00 

0.03 

0.68 

- 0.87 

166.10 

11530.16 15.26 

- 0.02 

226609 . 

0. 4014 

0.0000 

2626.305 

- 25.96 

65.16 

- 5.86 

0.50 

0.03 

11667.79 15.26 

- 0 . 90 

.. - 16623 . 

- 0.0 

- 0.1033 

26606760 . 

6191.7 

- 181.6 

16.3 

-1 .6 

187.3 

66.0 0.7 

- 16.0 

666 9 7 « 

657486 . 

6 . 0 

566732.76 

- 76.00 

- 0.00 

- 0.02 

- 0.23 

107.16 

0.65380 00 108.51 

0.90 

66. 

. 875.63 

12.01 . 

18.66 

- 3.69 

16.30 

9.03 

CG 

o 

sO 

_ 23 . 3 & 

8.06 0.00 

68. 12 

0 . 03 

68. 1 7 

0.0 












0 . 100 

6 . 5 

66.9 

r5.Q _ 

9.6 

_ _. - 76.2 . 

--- 0 . 6 . . 36912.1 

- 6.7 

- 638669 . . 

- 0. 71855 

0.50000 

0.0 

6077.8 

- 132. 7 

- 8.2 

26.6 

- 76.1 

0 . 6 . __ 1327.9 

.... -1.8 

102837 . 

. 3.0645 

... 0 . 0005 . 

90682.20 

90.00 

0.05 

0.36 

- 0.86 

166.21 

11536.96 15.26 

- .- 0.03 

225962 . 

0.3984 

0.0030 

2666.555 

- 23.71 

66 . 80 

- 8.09 

0.50 

0.07 

1 1329.65 _ 1 5 . 31 

- 1.75 

..... - 29092 . . 

. -0 . 0 

- 0. 1000 

26606760 . 

6182.6 

- 180.6 

2 3.3 

- 2.5 

188.2 

66.0 2.1 

- 13.9 

116703 . 

857505.* 

10.6 

566732.76 

- 76.00 

- 0.01 

- 0.36 

- 0.23 

186.07 

0.65660 00 108.53 

1.76 

1 89 . 

875. 63 

21.23 

18.63 

- 3.68 

16.86 

9.07 

75.92 

23.90 

8.06 0.00 

68 . 12 

0.00 

63.12 

0.0 












0.150 

6.7 

66.6 

.- 5 . 1 . 

9 . 1 

- 76.3 

0.6 95587.2 

- 6 . 8 

- 637277 . 

- 0.71726 

0.50330 

0.0 

6068. 6 

_ -10 3 . 1 _ 

r 7 . 6 _ 

20.5 

- 76,1 

... 0.6 . 3171.1 

- 2.6 

103096 . 

3.0411 

0. 0007 

90632.20 

90.01 

0 .06 

0.06 

- 0.77 

166.32 

11637.03 15.76 

-0.16 

776747- 

0 . 3063 

0.0000 

2662. 805 

- 21.68 

66.60 

- 7.15 

9.52 

. _ 0. 10 

11 19 6 , 32 . 15.61 

- 2 . 33 

- 41227 . 

- 0 . 0 

- Q . 1000 

26606760 . 

6173.6 

- 179.5 

19 . 2 

- 2.1 

1 83.5 

65.0 3.6 

- 13.8 

155426 . 

557530 . 

15.7 

566732. 76 

- 76.00 

- 0.05 

- 1.82 

- 0.23 

18 6 . 37 . 

0.655 10 00 10 8 . 72 . 

... 2.38 

689 ,. 

875.63 

30.45 


18.61 - 3.69 _ 14.96 9 -. 06 . 75.58 _ 23.94 8.06 0.00 _. 6 . 8.12 0.00 68.12 ___ Q.O 


0.200 

3.U 

44. 3 

-5. 3 

8. 5 

-76.4 

0.6 

1 80936.9 

-4.8 

-435603 . 

-n. 71 507 

0.50000 

0.0 

6059.5 

-183.5 

-6.7 

6.2 

-76.1 

0.6 

4058. 3 

. -2 . 5 

.103303. 

. .3.0177 

0.0039 

90482.20 

90. 01 

0.35 

-0*40 

-0.68 

166.42 

1 1684.38 

_.l 5 • 24. 

-0.24. 

. .224423 

. . 0.3923 

. 0. 0000. .. 

2481.055 

-19.27 

44.11 

-3.77 

0.54 

0. 13 

11117.72 

15.56 

-2.19 . 

-5 3321 . 

. -0.0 ... 

... -0.1000. 

24404760. 

6164. 7 

-178. 8 

4.9 

-0,6 

188. B 

46.0 

4.6 

-13.7 

214073. 

.. . 557554. 

.. 19.8 

566732.76 

-76.01 

-0. 2L 

-5-1.7 

-0.25 

184. 14 

0,65580 00 

108.93. . 

_ .2.40. 

.2463. 

_ 87 5.63... 

. 39.67 

_l 8*6.8 

-3.70 

LS..12 

2_.UA 

25*02 

24*3.0 

.3,£2 

0.00 


n. on 

63. 1 7 

0-0 


0.2 50 

11.2 

44.1 

- 5.5 

7.9 

- 76.5 

... Q ., 6 ,. 

221830.2 

0.0 

6050.3 

- 183 . 8 

- 5.6 

. 12 , J ... 

- 76.2 

. 0.5 

Al 5.7 . 8 . 

90482.20 

90.01 

3.02 

- 0.90 

- 0.57 

.166 , 52 

11725.62 

15,24 

2499. 305 

- 17.07 

44.01 

0 . 37 

0.56 

0.16 

11143.47 

15.74 

24404760 . 

6155.7 

- 178.9 

- 12.7 

1.3 

139.0 

46 . 0 

4 . 7 


566732 . 76 _ - 76 . 03 - 0 . 57 _ 

; . 1 8 . 55 - 3.72 15 .29 


- 9.07 


_^ 0. 28 1 . 8 . 6 , 26 . 0 . 6564 D 00 109.17 

JA , .49 , 2 - 6.06 8 • 07 _ 0 . 00 . 


- 4 , 2 _ 

- 433903 . 

- Q . 71 471 

0.50000 

- 1.6 

.. 103476 , 

. 2,994 7 . 

0. 0011 

- J 3 . 29 

223500 . 

... 0.3393 

0 . 0000 . 

- 0.97 

- 65781 . 

.- 0 . 0 ... 

- 0*1000 

J 1 L 3.7 

264482 . 

557565 . 

24.4 

1 . 54 .. 

.... 7007 . . 

. 875.63 

48 . 89 ! 

. 68 . 12 . 

0 .00 

6 3.12 

.0,0 


9.01 
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GOODYEAR AEROSPACE PROGRAM ZK836 
ASSESSMENT OF LOADS RESULTING FROM PAR AC MUTE DECELERATION SYSTFM 


CONVERTED FOR MSEC UNDER NASA CON TRACT NASd-29 1AA 

IHPUCI T'R£AL*8<A-H,0-Z) 1 " 

_ PEAL TM ( 302 ) « YA( 302 ) , >3(302 ) , YC ( 302) . Y~3 ( 302 ) , YE t 302 ) ,YF(302) t _ 

1 YGC302I ,Y8!302) .HEAOEmiAT 

DOUBLE' PRECISION L AMIDE ,L A.N20E , LA, MODE , NUDEG ,L AMO PD , KS , 

IK.SPKH2 , LTR, MJ, NUD/L HI , L H2, L i,C 2 , LOVl AN 1, LAN 2, LVMO,NU J_ 

‘ 2 , MUOE G , MUOO EG , N AP , NNP , L AM, L ANOE G , LS l ,LS 2 , L D P , L A MO P , N JP , N J P DE G 

3,M,M», tY, I YP.LT.L TD,L TC.NA ,NN,K,KSPK.Hl "3""' 

_A t MMA(l6>,MA t IlYP!l6),KBL,KBT,KHl,KH2,kRi,<KSjjj) 

0 I MENS l ONI A TMfiS ( 2 01» TJME IJ.il) . ARE A ! (20 ) 

_ COMMON' T , "of 7 x\ Z~, X P , I P , T HE , T IE P , XD , L 0 , XP 3 , Z P 0 , TH E D , THE P 3 , G A M , 3 A M P., 
_ 1ALP.ALPP, AM, AKP.DYPR, DY PKP , RHO , S , SP, D, DP , M, N?, [ Y, I YP, L T , L TO ,L TO , 
20CG,C,K,CA,CN,CM,CMQ,CCA{8 , lol ,CCN( 3,16) ,CC KH ,16) ,CCMQ(8 ,16) .CAP, 
_ 3CNP,CMP,CMQP, CCAP (8, 16) ,CCNP(B, 16),CCMP( 8, 16) .CCMwPl 8, 16) , V.VPjGRi 
AR, AA16.A) ,00(3, 3) , EE1 3) ,FF(3) .GX.OZ.OXP, OZP , QT HE , 0T HEP, APB AR , XtiAR , 

""S ZAAR, A AN (8 ) , A AMP I 8 ) , AAL PR I 16) . AAL PP E ( lfc) ,1 1 N . [OUT, “ _ 

60 ADT HE , 03DTHE , DAliTHP, OB Of HP, ASrtii , Bfl AR, A, 3, CHI , CHI 0, MU, MUD ,PH I , SP I i 

7P HID, PH I l ,PHI 2, BET 1.3ET2.EPS l, EPS?, KSPKH1 ,K SPKH2, LTR, OL , 

8S IGl.S IG2.ETA1, FTA2, TFNS.D AMP, S TH E , C THE » ST HEP V- THE P , LH 1 ,LH2 , _ 

“ 9A1BAR, A28A1 ,U1 BAR ,B2fl AR.BD 3 AN, AO BAK , L 0 ♦ L 1 , L 2, L AND , L AN l , L AN2 , NU, 3_ 
COMMON TOR, THEDU, THtOL.OTVC , T'.JRQ »BR I D » 1.4 M* T I ,CONF ,OTI ,CDAB J 

~ 1 . II u 1 61, .Si PJJI6i,lS 1 il S 2.LJP.LANOP.NUP.BPBAR.EP SP1, EPSP2,OLP.KS 

2,MMA,NA,TTG(«) , VVG18) ,VG,I I YP t KBL , KBJ ,KH l ,RH2 , T TENS.O ) ,K< S, _ 

_7"3CllNT(20),AALPtl6),AALPP(l6) _1__ 

1 FORMAT ( 3110.0) _ 

“ 2 FORMAT! 20A A ) 

50 FORMAT ( IHl, 9X, 'TIME ', 6X, 'RANGE • , 5X, •HORVEL* ,AX, • HORACC • ,AX, __ 

1* TUT4CC* , AX , * GAMMA • ,5 X, 'MACH N 0 ' 3 X , • T F N S l D N • , 3 X , *PHI • , 7X» ’OX* • _ 
2Tx7 i CA* f '8'X, *CAP*/lGx, 'TOROUE' , AX , • ALT l TUOE • ,2X ,' VERTVEL • ,‘3X , _ 

3* VERT ACC* ,3X, ' TOT ACCP ', 3X,' ■ GAN N AP * , AX,"' MAC H NOP v , IX , • 0 AMP ! NG* , _ 

A3X, 'PH l Q* ,6X, 'OZ' ,'dX, 'CN' » BX, • CNP' / 1 OX 7> MASS* , 6 X, I J H E TA • ,5 X , 

_ ' 5'THEVFL • , AX, • THFACC ', AX,* AX- S' ,6X, 'ALPHA' , 5X, 'OY-PR ' , 5X » ' LTO* 

67X , ' MU' ,8X»'QTHE' ,6X, 'CM* ,8X» ' C MP • / 10X, • N AS SP • , 5X» 'RANGcP ' , AX t 

7 • HORVE LPJ.,3 X J 1 QRACCP ' , 3X . « NCIR-G' , 5X, ' ALPHAP' , A X , * DY-P PP ' , A X , 

8 • LT • , 3X , ' M JD • « 7X , 'DXP ', 7X, 'CMQ« , 7X, 'CM3P • /l 0X,1 I Y* ,8X,.' ALU TUOE P' 

_ 91 X, 'VERTVEL P' ,2X , • VERT ACCP • ,2X, 8 AX-GP',5X, • VEL DC I T Y • , 2X , • DP • , 

_ 18X, 'LTD', 7X, 'CHI • , 7X, 'OZP* ,7X,' M,9X,'CAAP' /IDX ,1 1 Y_P' , 7X , • T HET API 

1 2AX, 'THFVELP* , 3X , • T HE ACCP • , 3X, 'NOR-DP *., AX, • VELOC I T YP * , IX, 'DEN SI TYI 

3,3X,'0CG',7X,'CHID',6X,'0THEP*,5X,'C',9X,'SP'/10X, 'LAMOA* 5X, 

A ' A ', 9X , »B « , 9X, 'PHI 1 « , 6X , ' P HI2» ,6X,'AiSAR' ,6X,«3BAR» ,6X,'NU P» ,7X. 

_ 5'APBAR' ,5X, 'BP BAR' , 5X , *LDP'»7X, 'GUST VEL'//) 

51 FORMAT! 8X.F8.3, 1X,8(F9. 1,1X1 ,F 9. 0 , lx, 2 ( F 8 v 5 ,2X )/_ 

18X ,9(F3 • 1 , 2 X ) , F 9 . 0 , 1X*2(F8.A, 2 X ) / 

28X,F3.2,2X, 8IF8.2 »2X) , F 1 0 . J , 2 I F§ .A,2J U 

“138X, F3. 3. 2X, 8! F 8. 2 1 2X) ,F 9.0, IX, 2! F8. 1, 2X)J 

A7 X.F9.0.2X, 8(F8.1,2X),F9.D, 1 X . F9 . J , IX , F 8 . 1/ 

59X , 61 F 9. 2 , 1 X ) ,E10.A,2 1F8.2 ,2X » . F9 .3 , I X , 2 ( FB .2 . 2X ) / 

68X, 12< F8.2, 2X )/ ) 

_52 FORMAT 1 ///2 OX, 'RUN ENDEO BY C ONSTRAINT'//) 

53 FORMAT(/iOX,8F10.3/10X,8Fl3/3_7) 

55 FORMAT ( / 10X , 'SHUTTLE _ROC<FL M3 T 3R » , 5X, 20A A) 

56 F0RMATI/10X,20AA) 
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5B FORMAT (IH , 10X , 'SPRING CONSTA NT ARRAY, <KS(8)»/) 

“59 FORMAT 1 1 H , 1 OX , 8 I F 1 0 . 1 ) / ) ' 

60 FORMAT ( 1H~', IPX, 'SPRING CON STANT TE NSION APKAY~, T TENS I 3 I « / > 

61 FORMAT 1 13A6.A2L I 

_ l I N = 5 _ _ J 

~ ' 1 O' IT* .6 

_ RE ADI r IN, 2 ) ATMOS 

RE AO ( I I N » 2 ) T I ME I 

READ! (IN, 2) 'ARE A l””"_ _ 

READ! I IN, 1) TT I , S~SP 1 .MBA, |~1Y P 

__ R E A 0 1 l I N , 1 ) T TO * V VO 

~ R E ADI I IN, 1 ) AAM , A AMP, A~ A L P F , A A L P P E_ " 

DO 20 1=1,16 

_ A At. P I I ) =AALPE ( I ) /5 7. 2 95 8 

20 AAL PP I I ) = A AL P PE 1 I )/57.2958 

__ READ! I I N, 1 ) I ICC A( 1 , J) , j = l , 16 ) , I = 1 ,3 I _ 

~ READ II IN, 1 ) I ( CCNI I , J ) , J = l , 16 ) , f=l ,8) _ 

R EAD! I I N ,_1 ) I ICCM I 1 , J ) , J = 1 , 16 ) , 1-1 , 8) 

READ! f IN, 1 mCCMO'l l , J )"," J = l , 16 » , I = 1,81 

“_P.EAOl I IN, l ) I ICC API I ,J ) , J=1 , iVllN’l, 8 1 

R'EADII IN,1 H ICCNPI I , J 1 , J=Jl', 16) i'l = 1,_8J 

READ! I IN, 1) I (CCMPI I , j ),j= l,ioi,I = f, 8 ) 

~ READ! I IN, 1 1 I (CCM JPi I, J ) , J = l, 16 ) , 1= l , 3J 

READ m N ,J ) TTENS , KK S 

105 PEAo'ii IN, 1) X,2, T HE , THF () , V ,GAM ,HH_H 

_ IFIX.FO.l.) GO TO 900 _ 

_J'“ READ! I IN, 1) THFP.GAMP ,VP,THEPD 

READ! I IN, if TOP , THEfjj , THE 0 lT» OJ VC_ 

READ! I IN, 1) APDAR.XBAk, 2BAK 

REA f) I I J_N ,_1 ) S,i),M, I Y,i TR,~C 

REAOUInTu DP.MP, lYP.DTl, Tl 

READ! I I N, 1) T ,!)T1 ,0T2 ,OTPl ,01H^2 , TOT C , T TT __ 

I READ! I IN, 1) L HI , LH2 , A IB AR , A23AR , 8 IBAR ,82 BAR 

READ I I I N, 1 ) LSi l LS2 1 DLP 1 <Hl L KHi 

READ! I IN, 61) HEADER _____ 

_RE ADI ] IN. 2 ) CONT 

THE=T HE/57 •2958 

GAM=GAM/57. 2958 ... ___ 

THED=TH ED/57.2 9 58 

THEP=T HEP/57.2958 

THEP0=THEPD/57.2958 

G A M P = G A M P/ 5 7 . 2 9 5 3 
XD=V*OCOS(GAM) 

2D=VTDS IN (CAM) ~ 

XPD=VP*DCOS (GAMP) 

2PD=VP*DSIN(GAMP) 

R = _ 6378377. I _ ' 

= 9.8 0591 60 








M A - 0 . 


CONF=0. 

11=0 


CALL BRIDLE 

WRI TE ! I OUT, 55) 

ATMOS 


WRITE! IOUT, 56) 

TIME! 



WRITE! I OUT f 53 J TTI : 

_WP I TE I IOUT, 56) AR E A I . 

WRITE! IOUT, 53) SSPI 

WRITE ( I OUT, 5(1) 

" WRITE! IOUT, 59) 1 KKS! II ,1=1,8) 

W RITE! IO-JT, AO) 

WR! TE (1 OUT , 59r!TTENS( I ) ,1 = 1 ,8) 

17 FORMAT ( 111 If 10 X, •AERUOYNAMIC P AR A ME TER S' / /) ■ 

18 FORMAT IlOX,' FOR BO OY ANGLE OF ATTA CK A RRA Y, AA LP!16) DEGREES* /> 

16 FORMAT! IX, 8F8. 3/) "" _ ~ ' _ " 2 J I 

3 FORMAT ( 2! IX , 8F8.3/J ) \ _"__22 

A FORMA T ( 1H 1 OX , * M ACH Nil MU Fjj AR K AY 2 A AM ( 8 » » / ) 

5 FORMAT! H , 1 OX , ' F URBOD Y AXIAL COEF . A RRAY , CCA tB.16)'/) 

“ 6 FORMAT 1 15 < 1 X , 8 ( F'i . 4 ) / ) j ' " ' 1 

7 FORMAT! 1H , 10X, • FORBODY NORMAL' COEF. ARRA Y, -C CN (~3, 16) '7 ) 

. 8 FORMAT! IH1, 10X, • FORBODY PI TCH MOM COEF. ARRAY,' CCM ! 8 , 1 6 ) •_/ ) _ 

9 FORMAT ! 1H , 10X, • FORBODY PITCH DAMPING COEF._ARRAY, CCM J ( 8, 16 ) »7) _ 

10 F OR M A T ( IH l x l 0 Xj '_A F T B 00 Y_A NGLE OF ATTACK ARRAY, A A L p P ( 1 6) DEGREE'/ 

11 FORMAT t IH ', iux, •AFTBODY MACH NUMBER AMR A Y, A AMP ( 8 ) •/ ) _ 

12 FORMATUH , 1 OX , ' AFT (TODY AXIAL COEF. ARRAY, CC AP ( 8 , 16 ) • /') ' 

13 FORMAT! IH , 10X, • AF T30DY NORMAL COEF. ARRAY, CCNP 13 , l 5 ) • / )_ 

2lA FORMAT i IH1 , 10X, 'AFTBODY P I TCH MOM ' COFF . ARRAY, CCMP( 8, 16) • 7) _ 

15 FORMATUH , 1 OX , • A FTBOOY P I TCH 0 AMP I N G COEF_. ARRAY, CCMQP(8, 16 ) « j ) 

_WR ITE! 1 OUT, 1 7) 

WRITE!! OUT ,19 ) 

WRITE! 1 OUT , 3 ) lAALPE(J) , j=l".,i_6) 

WRITE! IOJT, A) _ I ' 

WRI TE I l OUT ,16)1 A A M 1 J) , J = i , d ) 

WR ITE! IOUT, 5) _ _ 

iLE.11 F ! J.QU-Tt 6 me C A ( 1 1 J )., J = 1 , 1 6)7 1= 1 ,8).. — 

WRI TF ( IOUT, 7) _ 

WRITE! IOUT, 6) ! 1CCN! I, J).,J=I, 16) t I t-QJ 

WRITE! IOUT, B) _ 

WRI TE ( IOUT, 6) ( ( CCM! I , J) , J= l , 16 ) , I =1 , 8 ) _ 

WRITE! IOUT, 9) 

W RI TE ( I OUT, 6) ItCCMQU ,J), J=lTl6), 1=1,8) ’ ~ 

WR ITE ! IOUT, 10) __ _ 

WRITE! IOUT, 3) I AALPPE ( j ) , J= l, 16) ; 

WR I TE 1 1 OUT , 1 1 ) "" 22 ' ' 

WRITE! IOUT, 16) I A AMP! j'f,‘j = l,3) 

WRITE! IOUT, 12 ) " __ " 2 .” 2 - ^___ 

WRITE|J0UT,6J ! (CCAP! |",J( , J = 1 ,16 I ,t=1.8) 

WRITE! IOUT, 13) _ 

WRI TE ( I OUT, 6) (<CCNPil_» J),J=1 ,16), 1 = 1,8) 

WR ITE! IOUT, 1A) _ 

WRITE! IOUT. 6) IICCMPH.JI. J = 1 ,16). 1 = 1.8) 

WRITE ! IOUT, 15) 

WRITE! IOUT, 6) ( tCCMQP 1 I ,J) , J=1 , 16) ,1 = 1,8) • 





Page 59 
GER 15853 



_XPDD = FF(1| 

_ Z P 0 D= F F ( 2 ) 

NAP=( XPIjn*Cti<EP»ZPDO *Sf'H F P )/GR 

AIN P= 1 Z PP3 *C T H E P-X POD »STH | P \J_ GR 

_na= ( x 6 d*cthe+zoo^sthe )/GP ; ' 

NN= ( ZDU»C TH E-XDD»STHE )/GR 

THPOOO=FF ( 3 ) *j>7._2?58 

C A A P = CAP<‘SPI _ 

_VD=OSyRnXDD**2*ZDD**2J s 

_ _V P 0 = 0 S 0 H T ( X P D 0 * * 2 ♦ ZJP O D * * 2 ) 

W R I T£ ( I OUT , 51 ) r, X.XD, XPD. VO. GAMDE G . AM, T ENS , PH I DE 3 , QX. C A , C AP , TOR: 

IZ , ZD, ZOO, VPD.GAMPDE ,AMP,t)A *HP ,PH I ODE , Q Z , Ctt, CNP , M , T HEDEG, T HEDGE 
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2THE0DD.NA, ALPDEG, DY PR , L TO. HUDEG, 3 THE t CM, CMP , MP, XP, XPO , XPOD , NM , 

3ALPPJE ,DYPRP,LT , MUDDEG.QXP ,C -iOtCMOP, l Y,Z P, ZPD.ZPDO, NAP, V, i)P ,L TO, _ 
" 4CHIDEG, 3ZP, K,CAA<> , I YP.THPOEG.THPDOE , THPOOO , NNP « VP , RHO * DCG , CHIP OF , 
50THfPjC,SPI ,L AM DEG, A, H, PHt IDE, PHI 20 E \ AOAR, BB4R , NUPJ£G,APBAR,BP3AR, 
6L OP , VG - • - 

L = L ♦ 1 L_L__L______ 

T M < L r = r 

YA!L)=Z 2 1’ _ _ 

2 v B C L I = OYPR' " 1 _ 

YC ( L ) = V _ 

Y0( L ) = TENS l'_ “ 7 1 _ 

Y E I L ) = T HE DEG ^ 

_ YF(L)=THEDDE _ _ _ 

YG(L)=GAMOEG ' 2 ~ 

YH( L »= ALPDEG 2. “'LL LI LL. "L_ 

J J- 2 " " ' __ ■ " 

102 I F ( C JM S r J 103, 104, 200 

2 00 C ONTI MJ E 

CALI PL tk A JtTH.YA.YB, YC ,Y0 ,1,13,1, 4,2,5, HEAD ER) 

L CAL L PL Til AJ I TM, YE , YF . YG ,W.lVio7>, 1 1,6 , 8 .HE ADER) ~ 

GO TO 105 '2' 2 

104 WR I TE I I JUTjt'52'i." _ . 2 

CONST* l . 2 L 

j-n_ to i oi ' 

103 DO 74 J=l,4_ ; 

CALL SU OR * L 2. 

DO 75 1=1,3 

75 AAII.JI = EE! i)_*OT 

DO 76 1=4,5 _ 

76 AA I I , J J =FF(I-31»DT 


GO TO I 71, 72, 73, 741 jiJ 

71 X=X+XD*0T/2. . 

Z=Z+ZU*DT/2. ’ 2 

THE=THE*-rHE0*DT/2, ' 

XP = XP*-XP0*0T/2. 

Z P = Z°»ZP0<=D T/2. L 
T HE P = T HEP + f HFPD*DT/2. 
21_XD = X0»AA( 1,1) n. 

_ZD=ZDfAA!2, U/2. ’ 

_THED=THED«-AA(3, 1 )/2. 
XPD = XPi)+AA( 4,1) tZ. 

ZP D = Z P D»AA( 5, 1 )/2. 



J 


X P 0 = X P D— A A ( 4, 1)/2.»AA(4,2) /2. 

ZPO = ZPO'-AA(5,ll/2.+AA(5,2l/2_. 

~ T H E P D = T H E P 3 - A A ( 6 , l ) / 2 . #• A A ( 6 ,J ) /2. 

"_G0 TO 74 ’ ” 

) X = X»Or* (XD/2.-AAI1 ,l)/4.*AA(l ,2 f/4^) 

~ Z= Z*0r<‘ (Z0/2.-AA( 2, 1 > /4 . » A A ( 2,2 ) /4. ) 

THE = THE+0T*(THEI)/2.-AA(3,l)/4.*AAJ3»2 )/ 4. ) 

XP= XP-*OT*( XP0/2.-AA1 4 ,1 J/4.VaA(4,2)/4. ) 

ZP=ZP«-!H<M ZPO/2 .-AA< 5, n/4.+AA l 5, 2) /4. ) 

" THfcP=THEP*DT* I THtPO/2 .-AAI 6, l )/',l *-AA(h, 2 1 /4 . 
""XD = x:>-AA( 1, 2)/2. + AA(l,3) 

ZD=ZD- AA( 2, 2 ) / 2 . +A A12, 3 ) 

T H ED - Tit t: 0 - A A ( 3 j 2 ) / 2 • AM 3 , 3 ) _ _ 

2 XPO=XPU-AA(4,2)/2.*-AA(4,3> _ ~ 

' ZPO=ZPO-AA( 5, 2 ) / 2 . ♦ A A ( t> ,3 ) ”” ' ~ ] 

THE PO = THF P3-A A ( 6 , 2 ) / 2 • AA ( 6 ,3 ) ~ 

"T = T«H)T/2. II1I_ 

, 'CONTINUE ' ' ' 

XO = XO- A A ( 1,3) _ _ 

ZD=ZD-AA( 2, 3 ) 

THED=TMED-AA(3 * 3 J " 

XPQ=XP;)-AA( 4,3) _ 

ZPD=ZPU-AA( 5, 3) " _ _ 

THE PO -THE P U-A A ( 6,3 ) 

X= X— 5 T # ( XD + 4A ( 1,21/2. ) 

Z = Z-'JT*(ZD*-AA(2,2)/2.i 

THE = THE -DT* ( THEO*A4(3 ,2 )/ 2 . ) 

XP=X'’-l)T«(XPD + AA(4,2)/2.) 

Z P= ZP-D T* ( Z PD* AA( 5, 2. )/2.J> 

THE PMHFP-D T»( THE Pi) + AA( 6,2>/2. ) 

X = X +X D*')T *UAA (1,1) * A A (1,2) + AA1 1. 3) )«0T /6. 

_ Z = Z*-ZO*JTM A A ( 2 , 1 ) +AA (2,2) *;AA(2, 3 ) )*DT/6 . 

_ THE= THE *• THE i)*OT*< A A ( 3 , 1 )*A A( 3, 2) *AM 3.3 ) i »31 

XP = XP*-XPD*OT *• ( A A I 4 * 1 ) +AA(4,2)+AA(4, 3) ) T/6« 

_ Z P = Z P * Z PO * 0 T * ( A A ( 5 , 1 ) ♦ A A < 5 , 2 ) +AAI 5, 3 ) ) »p'r/ & . 
T HSP=THEP *-TH FP)*DT*-( AA( fa, j Lt M (6.2)*- A AI 6.31 ) 

XD = XD * ( AA ( 1 ,1 ) +2 . * ( A A ( 1 , 2 ) + A 4 ( 1 . 3 ) ) + AA ( 1 . 4 ) ] 

_ZD = ZJMAA(2 f l ) *• 2. * I A A (2 ,2) *AA(2, 3) ) + AAI2 ,4) I 
THED=THED+( AA( 3, l )*-2. * ( AA ( 3, 2 > *■ A A < 3, 3 ) ) » A4 ( l 
XPD = XPD*-( A A (4 ,1 )+2.»( AA(4 ,2 ) +A A ( 4, 3 ) ) 4- A A (4. i 
_ ZPO = ZPL>*-( A A ( 5, l )*2. *C AA ( 5, 2) *■ A A (JjjiJ )_*AA ( 5,4 
THE PO= T HE PD » ( A Alfa, 1 J »2 . *( A A [ 6 . 2 ) *-A A ( 6 . 3 ) ) *- 

IF ( T.LE.TQTC) GO TO 9? 

DT= DT2 _ 

DT P = DTP2 

GO . TO 99 . 

L CONTINUE.. 

CALL CALEND 

_ S TOP 

END_ , 

SUB ROUT INE SUBR __ 

_ IMPLICIT REAL»8(A-H, O -Z) 

_ DOUBLE PRECIS ION L AM 1 DE , L A M 2QE~, L A MODE , NUDE 3 , 
1KSPKH2 , LTR, MU ,MUO , LHl , LH2 , Ll,L2,LO,LAMl, L AM i 
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2, MU0EG,MJ0DEG,NAP,NNP,LAM f LAM BEG, L S 1.L S2.L0P.L A MOP, NUP , NU PDEG 

3, M.HP,i Y.IYP.LTtLTb.Lto'.NA.NN .K.K SV’KHl 

4, MMA( 16) ,MA ,1 I YP( I6j ,KBL,krtT,KHl ,KH2 , KPH IjKKS (J ) 

COMMON T, DT ,X,Z,XP, ZP, THE, THEP , XD , L l) , XPOj Z Pf) , THE 0 , THEPO , GAM, GAMP, _ 

1ALP, ALPP.AM.AMP.OYPR.OYPRPfRHj.S.SP.oJiPjMjMPLLY, 1 YP, L T , L TO , L T 0, 

2 n C G , C t < J c A , C N , C M .CMO , CCA( 0.16) . CCN1 8.16) .CCM18.16) ,CCM0 1 3 . 1 6 ) . C AP . 
3CNP.CMP, CM'JP.CCAIMa ,16) ,CCNP(8, lo) , C.CMPC 8 , 1 6) , CCM :)P 1 8 , 1 6 ) , V , VP , GR , 
4R , AA ( 6,4) ,30(3,3) ,EE(3) ,FF (3) , OX , OZ , OX P, OZP , QTHE , Q T HEP, APB AR , X 3 AR , 

5Z0AP, AAA( 8 ) , AA MP( 8 ), AAL PE ( 16) , AALPPEt 16) , i IN, IOUT, 

60 AD THE , DRDTHE , 0 AOTHP, DfiOT HP, A8 AR , OH AR, A, B, CH I , CH I 0, M J, M JD , PH I , S? I , 
7PH[0,PHI l, PH I 2, BE T l,3ET2,EPSl , EPS2 , KSI’KHI ,KSPKH2 , LTR ,DL, 

8 S I G l , SJ G2 , ETA 1 , FT A2, T ENS t IT AMP, S THE , C THF , STHEP , : THE P , LH 1 ,LH2 , 

9 A IB AR , A2 i) AR , P 1 OAR , 82 BAR ,U08 AR , A3 BAR , L 3 , L l . L 2 , L AMO , L AM 1 , L A.M2 , NL) , G_ 
COMMON TOR, THFJU, TilFOL, OTVC ,TORQ ,BR [ 0 , LAM, Tl ,CONF ,OTI ,COAU """ 

1 , TT I ( 16 ) , SSP1 ( 16 ) ,LS1 ,IS2, LOP, l AMOP,NUP, BP BAR , EP SP 1 , FP SP2 , 3LP , K S 

2 , MMA , M A , TT > l 8 ) , VVG( 8) , VG, I I YP , KBL , K BT , KH l , K H2 , T T ENS ( 3 ) , KKS , _ 
3C0NTI 20), AALPI 16), AALPP116) 

J I I =1 

G = GR*(R/(Z«-R) I **2 


” RE AC T I ON CONTROL SYSTEM' 

IF (TORO. EG. 0. ) GO TO 49 ~ 

IF ( T HE 0 Ai f . TH F QL) GO TO 48 

I F ( THEO.LT. -THEOL) GO TO 47 

IF (T-TC.GT .OTVC) GO TO 40 _ _______ 

IF ( THEO.GT.O. ) TORQ=-TOR* { i •- ( T-TC ) /JDT VCJ 

IFITHED.LT. U. ) T3RQ=T0K*( i.-( TrTCL/DIVC) 

GO TO 42 _ 

48 TOR 0=^ TOR _____ 

TC=T ~ 

GO TO 42 

_47 TORQ = TOR _ 

TC= T ~~ I 

GO TO 42 . 

49 I F ( TH EO.LT. THFOU.AND. THEO. GT.-THEOO) GO TO 40 

_ . I f7~T HEO.GT .THEO U) GO TO _41 _ 

TOKU=TOR 

TC=0. 

_ GO TO 42 _ 

41 TORQ=- T OR 

_ TC = 0. 

GO TO 42 

40 TORQ = 0. 

42 CONTINUE 

T ORQ=0 .0 

1=2 

400 IF ( T. LE .TTGI I ) ) GO TO 401 

l=IU 

GO TO 400 

401. TSL=( T-TT.GJ.I- l ) )/( TTGI 1 )- T TG ( I - 1) ) 

VG=VVG( l-l ) +(VVG( 1 i-VVG(I -l) )*TSL 

_ CALL OENSM( Z, PR, R HO, VS) 

C CALL DENS (Z.PR.RHQ.VS ) 
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V = 0 SOR T(XD* *2* ZD ** 2, _ 

v p= Dsg'a t Yx Pb** 2 «- 2 Po** 2 ) 

~A M =0 S Oil r ( ( X 0- VG ) * *2 ♦ Z 6**2 » /vs 

AMP = DSO^T(l XPO-VG)**2«-ZP3**2) /Vi 

DYPR= .5*RH0*< UD-VGI?*2 +20**2) 

QYFKP=. 5«R-I0»( tX Pl)-VG)»»P»2Pr)*»? I 

S r HF= ns IN ( T HE ) 

C TH E= DCOS ( T HE ) I " 

' STHEP= 1 SIN( THEP) 

^CTME-’OCOSI TtlFP ) ' 

GAM = I1ATAN2 ( 20, XO) " ' * 

C A MP= D A T A N 2 ( Z P[),XPI ) j 

ANG= P AT AN? ( ZD, (XD-VGI)_ ___ 

I F ( ANG .LT.O.O) A.'IG* A NG + 5.2331854 

201 THEA=6.233ld6+THfc J 

_ A L P = T HE A- ANG 

IF (ALP.GT.3.1 ‘.1592 71 A L P = AL P -6 . 28 3 135 ~ 

A MGP- (JATAV212PI) , ( X°0-VG) ) 

__ I F ( AN G P . L T . U 1 0 V A N GP = AN GP ♦ 6 - 23 3 1 8 54 

'3_0i”THEPA = 6. 283185*- THEP ' 

ALPP=ri)EPA-ANGP ‘ 

IF ( ALPP.GT. 3. 1415927J ALPP= ALPP-6 .2 83185 
700 CALL AERO __ _ 


C SPIOLF, RISER, SUSPENS ION GEOMETRY 

C_ _ 

A = XP + APBAR* CT HE P-BP BAREST ME P-X -AR AR*C TH E »B BAR» STH £ 

8= ZP +APB AR*STHEP+aP0AR*CI3EP-Z-ARAR*ST HE -BlLAR*C THE 

A D= X o 0- AP B AR * T I i E P 0 «STH EP-BP B A~R » T H E PO*C TO E ~P - X3 + A 3 AR* THEP* S T HE + 

1BB AR» TH EP*CTIIE 

B 0 = Z P 0 + APB AR*TMEPO*C THE P^-B 11 B AR* THE£2* STHE P-ZD -AB A R»T HED*CTHE» 

1BBAR*THE0*STHE 

1' lT = 0S0!<T( a**2 + 3**2’] _ 

IFILT.LT.LTO) GO TO 35, 

TENS=K* (LT-LTQ) _ 

GO TO 36 

3.5 _J ENS = 0 . , 

_ _36 CH I =0 A T AN ( A/ B 1 

^ MU=-1.570 79 63-TH E P-NU P-CHI 

PHI = 1 .570 796 3-T HE-L AM- CHI ~ 

IFIB.LT.O.) MU- 1 . 570796 3— THE P-NU P-C H I 

IFIB .LT.O.) PH 1 = 4. 7 123 88 9-THE -LAM-CHI 

20 PH i B=PH I 

‘ i =2 ' _ ~ _ ' . 

600 I F ( TFNS. LE.TTE NSt I ) j GO TO 601 

I = 1+1 

GO TO 600 

601 T E NS L = I TENS-TTENS 1 1-lTl/I T TENS ( I l-TTENSI I - 1 > ) 

K S = KKS ( 1-1 ) ♦ ( KKS ( I 1- KK S 11- 1 ) ) * T ENSL 

K SPKH 1= ( 2. *KS*KH l j / 1 2. * KS+KH 1 ) 

KSPKH2= 12 .*KS*KH2 ) / 12 ,»K S+KH2 ) 

IF (PHI .GT. PHI 2 • OR. PHI .IT. -PH 1 1) GO TO 15 

PL 1= ( TENS/(XH1«0SIN( SIG1»S1G2) ) )»"(-DSINIP-II ) » OCOS t S I G2 ) 

1 + DCOS (PHI) *DS IN ( S IG2 ) ) 
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t>L 2= ( TENS/(KH2*r)SlN< SIG 1 + SIG2) ) ) »(DSI N( PHI ) »DC 0S (S IG1 ) 

1~'»0C0S( PHI ) *DSINl S 1G1 11 

LHI = LH1»3L1 __ I 

LH2 = LH2«-DL2 

CALL BRIDLE _ 

LHl = LH l- ) L 1 

_ IM2 = LH2-0L2 _ 

’ ORA R = 6 J JAR _ " _ 

A04R=A0BAR ~ ' _2_2_ ~ _ _ 

PH I L = S I G1 2_1 21_ I 1_ 

PH 1 2 = S I G2 ' . _2~ ~ ■ _ 

KP H I = K 3 L * S BT/<KBT»(r>SOS(PHI )T« *2 » KB L » ( PS I N I PHI ) )*»2) 

K=2.*KS*KPHI/ (2.*KS+KPHI) _ 

lto=ltr __ __ 

_ LAM=NU _ _ _ I 

GO TO U ” " ’ 

I 15 mPHl.LT.-PHIl ) GO T 0 lb _1 I _ 

A B A R = A23AR 

BBAR=B2 BAR _ ___ 

” L TO = L TR »LH2 

K = KS°KH2 

_ PH12=ETA2 ' 22" _ ' 

LAM=L AM 2 ' ~ 

I £ ( LA 0. 01.3.14150271 LAM^ L AM 2- 6. 2H3 1 8531 

_____ GO T 3 14 __ 

16 A BAR= A1 BAR _ 2_ 1 ~ 2 " _ ~ ' 

I 222 1 B B A R = B1 BAR 21 

LT0 = LTR*-LH1 _ _ 

K = KSPKHl ~ ' 

_______ PH I 1 -5 1 A l 1 

L4M = L AM 1 _ 

14 OLS l = ( T ENS/ I KS*OS I N ( EPS PI +EPSP2 J ) I «( DC OS ( EP SP2) »DCO SI M U»MUP ) - 

IDS I Ml EPSP2I *0SIN( MU*N1JP) j __2 ' 

DL S2 = ITENS/(KS»0SIN(EPSPj» EPS P2) ) 1 * ( OC OS ( E P SP 1 1 * PCD S 1 MU»MU P ) + 

10SINIEPSP1 )*DSIN(MU*NUP)I 

_L S 1 = L 5 IOLS 1 

LS2 = LS2«-n>LS2 

C I. CALL SUSPEN _2 2 ' 

_ LS1 = LSI-3LS1 ; 

LS2 = LS 2-OLS 2 _ _ 

_ A=XP* A P3AR*CT HEP-BP BAR* STHEP-X-ABAR *CTHE*BBAR*S THE 

B = Z P»APBA3»STHEP»BP3AR»CfHEP-Z-AHAR*STH5-BBAP»CTHE~ 

_ AD = XPr)-APBAR*THEPP*STHEP -OP BA R»THEPD*CTHEP-XO»AaAR»TrIED»STHE» 

1BBAR*1HFD*: THE ’ ”2 ' 

B 0= 2P 0 * AP BA R * T HEP 0*0 TH E P- B P B A R * TH E P0» SHE P- 20- A 3 A R * T HE 0*C TME* 

1BBAR*THED*STHE 2 1 ~~~ 

LT=OSQRT( A»»2»B**2) 2 2 

L T 0= 1 A* AD+B * BD 1 /L T - 

CHI =0 AT AN ( A /B 1 

■ MU=- 1 . 5 707963-THEP-NUP-CHI 

2 PHI = 1 .5707963 -THE - L AM - CH I 

IFIB.LT.O.) MU- 1 . 5707963-TH EP-NUP-CHI 

IF(B.LT.0.)PHI = 4. 7 123 8 3 9- T HE -LAM-CHI 

CHI D= [ B*AD-A*BU)/ (LT»*2 } 
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MUD=-THEPD-CH 10 

PH I 0=-I HE 0- CHID _ I_ ~ 

' IF ( IT. L r.LTO) GO T O 30 

TENS = K * I L T - L TO J I 1 

DAMP=C*LTD 

GO T O 31 ' ' 

T EN S= 0. __ 

D A M P = 0 • _ 

1 1 1 = 1 1 in 

I F l III .GT.i 0 )" WR I TF.I I OUT, 50) 

IF( III. Of. 10) GO TO 19 

_f n '■) M A r ( 2 0 X , 1 1 TCRAT10N DOES NOT CONVERGE*) 

JFI UAHS (DHI-PHIB) .GT. . 0 083 ) _GJ _T3_20 

DA3T.Hc=AllAri*S THE *• BBAR*C THE _ ~ 

000 rilE=-AHAR*C T H E ♦ rl U A R * S7HE ' I 

0 A 0 T H P a - A P 0 AR * S T H f.V -U P B AR * C t H E P 

0 BO Til P = APB AR* CTHf P-BPBAR*S THFP " 

_C AL L _M A TR_I X " " 

FF (!) = (-( TCNSfDAMP ) *A/L T + CXP I/HP 


I_f FF I 2 I = I -( TE NShOAMP )*3/LTMJZ Pj/ MP-G^MP -M A ) /MP L_ I 

_~I FF ( 3 ) = (-{ TFNS+OA'IP )*( A*DA3 T'1P» B *PB3 THP I / L T «■ QT HE P ) / I YP ’ 

1 '^RETURN ~ I_I 

II 'END I"‘ ' . "I”"_II _ 

S U BROU TINE AERO 

IMPLICIT RFAL*HIA-H,0-ZI 

OO'JRLE PR EC I SI UN I AMIDE , L A M2DE , L AiMOOE ,N UOEG , L A~MO PO , KS , 

„ l K S P KH2 , L T R , M J , M J 0 , L HI, L H2 , QIIT , L 0 1 L.AM1 »JL AM 2 1 L A MO,~N U" I I 

2 , MUI1E G , MUD3 EG , N A P , NNP , L AM , L AMOE G , LS I , l S 2 , U P . L A MO P , NJPj,NJJPDEG 

3 , M , M ° , I Y , I YP , L T , 1. TO ,L TO ,NA , NNj_K » KSPKH 1 I 

4, MM A I 16 ) , MA, I IYP( 16 ) , K f j L , K 0 T , KH1 , KH2.KPH I , < K S I B ) 

COMMON T » 3 T , X , Z » X P , ZP , T HE- , THEP , XO , ZO , X PD , ZP Q, T HEP, THEPO, GA M, GAMP, 

JALP, ALP 0 , AM, AMP, UVPR, DYPRP S, SP.,q.,O.P ,MjMP_, I Y , I YP, LT,tTO .LT D. 

20CG ,C , K ,C A , C.N,C.M, CMQ, C C A C 8 ,16) , CCN ( 8 , 16 ) , CCM t 8, 1 6 ) , CCM J ( ft , 16 ) .CAP. 

3C NP ,CMP , C M3 P , CC AP t U , 1 6) ,CC NP ( 8 , 1 6) , CCHPI 8 ,1 6 ) ,CCMO° (3 , 16 ) , V , V P, GR, 

4R , A A ( 6 , 4 ) » Dpi 3, 3), EE I 3 ) , EF ( 3 ) , (J X , QZ , QXP , QZ ? ,3 THE , 3THE P , A PO AR, XBAR,_ 

5ZJ3AP , A AMI R ) , AAMP( 8 ) , A AL PE ( 16 ) . AALPPE 11 6 ) . II N. IOJT . 

6DADT HE , 000 THE , D AO THP , DQOTHP , AB AR ,00 AP, ,A , B , CHI , CHID,. MU, MUD, PH I jSPI 

7 P H I 0 , P H 1 1 , P H I 2 , OCT 1 ,BE_T2,FPS1, EPS ~2, <SP <H 1 . < SPKH2. L TR, DL . _ _ 

8S I G l , S I G2 , c TA 1 , E T A 2 , TE NS , 0 AM P , S THE , C THE ,.S f HE P , CTH EpYlHI , i H2 ,„ “ 

9A 1 BAR , A2BAR , B 1 i3AR , 0 20 AR , 000 AR , AOBAR ,L 0 , L 1 , L 2 ,L AMO, L AMI , LAM2, NU,G 

COMMON TOR, THEOU,THEOl,UTVC,TORg, BRIO, LAM, TIjCONF, DTI, COAB 

1. TT II 16). SSPII 16) .LSI li S2.L0P.LAM0P.NUP. BPBAR .EPS PI . EPSP? .01 P. KS 

2,MMA,MA,TTGI8I,VYGIR),VG, I 1 Y P , K8L . KBT , K H 1. < H2. T TEN S I 8) .KK S. 

3CC1NTI2JI , AALP 1 1 6 1 .AALPPI 16 ) 

1 = 2 _ " I_ ' 

5 00 I F I A M . L E .A AM i I_£l_ GO fo 50_1 I 

I =1 > L .'I 

G 0 TO 5 0 0 

501 A M S L = I AM- A AM I L; U_]_/ I A AM 1 1 ) - A AM (1-11) 

J =2 I 

60 0 1 F I PASS I ALP ) . L E.AALPI J ) ) GO T3 601 

_ J = J + 1 

GO TO _600 

601 AL PSL = I DABS I AL P l-AALP I J— 1 ) )/ I A AL P t J I -AAL P I J - 1 ) ) 


Page 66 
GER 15853 



C A= ( CC A ( I, J-l )-CCA( I-_l, J^.1) >*AMSL*CCA(I-1 , Hit ( tCCAJI , J)-_ 
lCCAt I-l , J ) ) *AMSL +CCAI I-l, j )-< ( CC At I , J-U-CC At t-1, J-l I )*AMSL«- 
2CCAI t-l,J-l)> )*ALPSL L J 7 ~~ ~~ 

cn= ( c c n! i i, j-i >-ccni i - i,j-i ) )’*amsl«-ccmi ir_i» Jti>«-i tccNti , j>- 

lCCNt I -1 , J) > *AMSL*CCN( 1-1, J )-( (CCN t I , j-l 1 -CCN ( l_-l, J-l ) > ■S'A M SL ♦ 
_2C C N < t - 1 , J - 1 ) ) > * A L Pit . 

c«=tcC'iu,j-n-cr.M( i-i, j-_ij )*amsl ♦ccmij- i,j-h*( (ccmi t , j j - 

ICC M I I -l ,J) ) *AMSl»CCM( I — l»J I — t tCCMC I , J~-l ) -CC.M ( l -1 , J-l ) ) -AMSL ♦ 
2CCT t i-i, j-i m *al.p sl 


CM'JrtCCIJl I , J-l J-CCMOt I-l , J-l 1 ) * A .MS L ♦ C C M 
JCCMOt I-l, JJ ) ^A'lSL + C C M 0 { I - 1 , J ) - t ( C C M Q ( I , J 

2CC MO (7 -1 * J-l" TT ) *A L PS L _ __ 

IFIALP.LT. O.) CN=-CN ” _ 

IFtALP.LT.O.) CM= —CM ' 

1=2 _ ' j 

) !F( AMP.LE.AAMPdJ > GO TO 701 

_ l = I ♦ L 

GO TO 700 _ 

L AM PSL= t AM P-A AM°( I-l ) ) /( A AM P ( [j-AAM P( I- l )_ 

_ J=2 _ 

3 I F ( DAO S( A LPP ) . LE . AA LP P j Jl ) GO T~0 80 1 

J = J + l _ 

GO TO 800 


L ALPPSL=(OABS( ALPP)-AALPP(J-l) )/(AALPPt J) 
CAP=ICCAPl I ,J-l l-CCAP ( I-i,J-ij )*AMP$L +CC 
ICC API I-l, J) )*AMPSL*CC APlI-l,J|-( ICC API I, 
2CCAPM-1, J- 1) ) I^ALPPSL 
CNP=(CCNP< I , J-l J-CCNP < I-l, J- j) I *AM~PSL«-CC 
1C C N P t J - l ,_J ) ) * MPS L + CSNPtl-l, j)-( ICCNPt I, 
2CCNPI I-I, j- H 1 ) *ALPPSL 
CMP=(CC MPtI , J-l I -CCMPI I-l , J-l ) ) *AMPSL+CC 
lCCMPt I-l, J) )*AMPSL«-C:NPII-i,j»-i iCCMPII, 

2CCMP ( I-l, J-l) ) ) *Al P PS L ___ _ 

C MP=r, MP+ APB AR4-CNP/0P _ __ 

__CM Q P= (CCMOPiLi Jjr. I I - c: K QP t I - 1 . j- 1 ) ) * AMPSL 
It tCCMOPt I, j )-CCMOP ( I-l, J) ) *AMPSL KCMOPt I 
2CCMQP( I-l , J-l ) ) * A M P S L +CCMQPI I.-ljJ r l ) J J*A 

_ I F ( ALPP.LT.O. )CNP=-CMP 

IF( ALPP.LT.O. )CMP=-CMP 

QX = -OYPR*SMCA*CTHt;«-CN*STHEt 

0 Z^D Y PR* S * i.CN *Q THF-CA»STHE) 

QTH£=JYPR*S*0* ( CM+CMQ*T HE 0*0/ VJ +IQRQ 

1 = 2 

3 IFtT.LE.TTltm GO TO 901 

1 = 1+1 _ _ __ 

GO TO 900 _ 

i TJ _S_L=( T-T TH 1-1) )/(TT 11 I)-TTI( I-l)) 

_ SP l = SSPIt I-l»ti SSPIt I l-SSPI (I- l) )*TISL 

M P= M P -MA 

_ MA = MM A ( I-l ) + (MMA( I )-MMA ( I- 1 ) ) *Ti_sX 

_lYP = IIYP(I-l) + tIIYPtI ) - 1 1 YP t I-l I )*T ISL 

MP= MP FM A__ ~1_ Z I j 

IFlT.LT.ofl) CNP = CNP *T / DT I 


Qtl-l, J-l )♦( (CCMO 
-1) -CCMtJt I-l , J-l ) 


t l,J )- 
) »AMSL+ 


-AALPPtJ.nl) I 

API I-l, J-i) 4-X C CC A F 
J-l ) - C CAP t I.-l , J-l ! 

nT P 77 - 1 ,7 - 1 > 7 ( ( c C M F 

J-l )-CCNP( 1-1, J-l I 

MPU-l, j-l)7t (C CM f 
j-i )- cc mp i i -i , j-i i 


FCCMQPt I-l , J-l )+■ 
-1, J >-( t CCMQPt I , J- 
LPPSL 
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c 


_ 1F(T.LT. DTI) CMP=CMP*T/OTl 

IFtT.LT.DTIl C‘TLP=CMji>»f/D~T 1 ' 

QXP=-OYPRP*(SPI *CAP*CTHEP*S~P*CNP*STHE_PJ 

OZP = DYPpP*tSPI *CN P *C T MF P- S.PI *CAP*ST HF P ) 

QThEP=OYPRP*S PI »UP* tCMP»CMUP«rHEPO»DP /VP 1 

rf turn 

_ END _ 

_ SUBROUTINE MATRIX 

1 IMPLICIT REAL *3( A-H,3-Z ) 

003 BL F PRECISION L AM 1 PE ,L A M2QE . L A MODE . N J DEG . L AMOPO. KS. 

IKSPKH2.LTR, MU.MUJ ,LHl , LH2 f L 1 ,_L 2 j.L.0 j. L A Ml , LA M2 ,_L A MO, N U 

2,MUPEG.HJD0EG,NAP. NN n . L AM, LAMDEG.LS1.LS2.L0P.LAM0P.NUP.N1JP0E3 

3.M, MP, I Y, I YP, LT rb.LTO ,NA,NN, K.KSPKH1 

A.M'-IAt 16 1 ,MA , I I YPI 16) .KBL.KBT ,KHl ,KH2 ,KPH 1 ,KKS ( 3 ) _ 

COMMON T.DT.X.Z,. XP.ZP.ThL, THEP , X0 , Z D , XPD , ZPD , THED , THEPO ,GAM , G AMP, 
1ALP, ALPP, AM , AMP, UYPK, OYPRP , RHO , S , S P , D , DP , M , MP , j Y , I YP, LT, LTD, LTD, 
2PCG, C.K.CA, CN,CM,CMO, CCAI 8 , 16) ,CCN( 3 ,16) ,CCM( d, 16) ,CCMQ(8 ,16) ,CA», 

3CNP , C M P , C M;jp , C C APJ ;! , 1 6.) . CC N P ( 3 . 1 6 ) . CC M P t 8 , 16) , C C M 0 1> ( B . 161 « V . VP . G.R . 

AR , A A ( 6,6) , 3UI 3*3) , ELI 3) ,FF ( 3) , QX , QZ , OXP, QZP , OTHL , QT HEP, APB AR , X3AR , 

_ . 5ZBAR, A \Mt 8 ) ,AAMP( 8 ) , A ALPFI 16) , AALPPEt 16) I IN, IOUT, __ __ 

60 AD THE .DBOTHF .OADTHP, DUOT H P, A3 AR , BBAK , A , B , C H I , CM I D, MJ, MOO, PHI , SP I , 

7PHID.PH II, PIII2.QET1, BET2.EPS l,FPS2,KSPKNl,KSPKH 2, LTR,DL, 

8S IGl.S IG2, ETA1, ET A2.T FNS, DAMP, SThF , C I HE , STHCP , Z THE P , LH 1 , LH2 , 

9AJJJA R , A2 BAR. , 81 BAR , 8 2 BAR ,!) 03AR, AO BAR , L0.L1 , L2 . L AMJ . LAM 1 , L AM 2 , NU.G 

COMMON TOR, THFOU, THEDL , l~)f VC , T U «0 , 3R I D , LA M , T I , CONE , DT t , CDAB 

1 , TT II 16 ) , SSPI (16 ) ,LS1 ,LS2, LOP, LAMOP, NUP, BP BAR, E P S P 1 , E P S P 2j 5 L P , K 5 

2, MMA,MA,TTj(8) ,VVG(8I ,VG,I I YP ,.KBL tKBTj_KHi ,KH2 ,JT E_Ni(3) ,KKS , 

3CONT(20),AALP( lfcl, AALPP.U6)__ 

00(1,1 ) =M 

00 ( J , 2 ) =0. 

001 1,3 )=-M»I ZBAR*CTHS +XBAR«STH E 1 

not 2 , 1 ) =o, 

DD t 2, 2 ).=M 

0012,3) = -M* ( Z BAR* STHE-X BA R *CTHE 1 

00(3,11=00(1,3) 

DP(3, 21=00 (2.31 

00 ( 3 , 3 ) = ( Y 

THE 02 = THE 0**2 

F F ( 11 = 00(2, 3)*THE02*( TENS + P AMP ) »A / L T» 0 X 

EE ( 2 ) =-00 ( 1 , 3 ) * THE 02 ♦ ( TENS *0 AMP ) »B /LT- M*G+QZ 

E E ( 3 ) = -G*DJ ( 2 *3 ) - < E ENS»DA MP)»(A*DADTHE+B*0BDTHE)/LT+QTHE 

E PS = . 1 o-ll ; 

CALL CR0UT(D0,EE,3 ,3, E PS, I ER SW ) 

I F ( IERSW.EO-0) RETURN 

WR I TE ( I OUT , 50 ) 

STOP 

j> OF 0 P M A T( »// 2 0 X jJLl NCUNS1STENT EQUATIONS*) ; ; 

END. 

SUBROUTINE BR I OLE 

l M P L I C I T R E AL* 3( A-H, O-ZJ 

DOUBLE PRECIS ION L AMI DE ,L AM2DE , L AMODE , NL) DEG . L AMOP D . KS . 

IK SP KH 2, LTR , MU , MUO,LH 1 ,LH2 , LI , L 2 , LO , L AMI , LAM2 . L AMO , NU 

2, MUD EG, MU 00 EG , NAP , NNP , LAM, LAilPtC, LS 1 . L S2 . LOP . L A MOP . NUP .NUPDES 

3 , M , MP , I Y , I YP , L T , L TO , L TO ,NA ,NN , K , KS PKH1 


O ^ 
M &> 
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4,MMA<16),MA,I IYPU6), K BL , K fU , < KH2 , KPJ[U KKSJ 8 ) 

COMMON T, OT , X , Z , X P . ZP , T EE ,TH E P , X D, ZD ,X PD .ZP O.THE D,_T HE P D , GAM , GAM P , __ 

1 AL P* AL PP» AM, AMP,DYPR,DYPRP , RH J , S , SP ,[) ,DP ,M, MP,l Y, I_YP,LT,L TO ,LTO,_ 

2DCG.C, K.CA.CN ,CM ,CMU, CCAI 9 , 16 ) , <7CN ( 8 , 16 ) , CCM ( 8, 1 6 ) , CCM3 1 B * 1 6 ) , C AP, 

3CNP,CMP ,;M3P,CCAP( R, 16) ,CCNP(8, 16) ,CCMPJ 8,1 6 ) t CCMOP (3 ,16) , V , Y P » GR , 

69, A A ( 6 , 6 ) t OO( 3, 3 ) , E E ( 3) ,F F ( 3) . QxV.l Z . 0 XP . (JZP . 0 THE , □ THE P . APB AR . X :i AR . " 

5ZRAR»AAM(8) , A AMP I 3 ) , A AL PEI 16 ) , AAL PPE ( 16 ) ,1 I N, I OUT , 

~ 60A0THE , 03 D THE , DAO TUP, OBDTHP, AJAR ,8 BAR , A, 6, CM I , C H I 0 , MU , MU D , PH I , S P I , 

7PH1D, PHI 1 , PHI 2, DEI 1,8 ET 2, EPS I, EPS 2, KSPKH 1 , K SPKH 2, L TR, DL, _ 

_ 8S I G 1 * S I G2 * E TA1 » £ T A2 » T ENS » DAMP , S THE , C T HE , $ T HEP , CTHEP , LI II , L H2 , 

9A 1 BAR , A2RA1 , B l.iAl .I12BAR .BOHAR , A0R1G ,LO,L 1 ,L?,LAMO,LAMl ,LAM2 ,NU,G 

COMMON TOR , THE UU.THE PL, PTY C , TU *0 . OR ID. L AM, T 1 , CONF , D J I , CD AB 

1, TT I ( 16) , SSPI I 1*6) ,LS1 ,L S 2 , LOP , L AMOP , NUP, BPB AR , EPSP1 , EPS P2 , f)LP , KS 

2, MM A, MA,TrG(ft ) , VYG I 9) ,VG, l l YP, <BL .KBT.KHl ,KH2 , TTENSl 9) ,KKS, 

_ 3C0NT(20),AALPI16) .AALPP116) _ 

'53 FORMAT! 1H1, 15X,20A6///I _ _ 

56 FORMAT C 1 H ,10X, ‘INITIAL YALUES, ENGINEERING UNITS ARE METRIC 1 METE 

IRS , N-: W TON , SEC ) •/) 

5 5 FORMAT ( IH ,/lOX, * L H 1 ' ,7X, * L H 2 ' , 7X , ■ A 1 BAR * , 5 X , • A2BAR* , 5X, • B1 BAR • ,5X , 

1 ' B2 BAR * , 5X « • LI ' , BX , ' L 2 • ,3 X , ' OL ' , HX , *L0',8X , • BE T l • , 6X, 'BE T 2 • /10X, 

2'EPSl* ,6X, *EPS2* , 6X, 'ETA1* ,6X, • ETA2' t 6X, JSIGl* ,6X SIG2 • ,6X , _ __ 

3‘LAML • , 5X , ‘LAMP* , 6X, ' L AMD* , 6X, • N U • , 8X, • A 03 A R • , 5 X , • 8 03 AR • / 1 OX , 

6 ' K S PK8 1 ,5X»' KSPKH 1 ■ , 6 X , • K S PKH2 ' , 6X , • T- I N F ' « 5X , • DT - l NF • , 6X » _ _ 

5 * TH E DU «j 5 X, • J H E DLJa.5X.tl TORQ UE' . 6X OT- VALVE « . 2 X.»LS1» . 7 _X . ‘LS 2» . 

_ 67 X, • DL P* / 1 OX, ' LOP ' , 7X t 1 LAMOP • 1 5X, 'NUPJ , 7X,_?AP-3AR • ,_5X, ' BPB A R • , 5X, 

_____ 7* E P S IP ' , 5X, ' E PS2P • , 5X , ' DT‘ , 8 X , • KHl 1,J7X , • KH2 ' / )_ ’ 

56 FORM AT (3X, 1 2 ! F8 . 3 , 2X ) /8 X, I2(F3.3,2X) / 8 X t 3 ( F 9. 0 » 1 X) ,9<F8.3,2X) 

1/8X,8(F3. 3,2X1 ,6<F9.0jlX)// ) _ 

C __ 

C STATEMENT FUNCTION 

D ARCOS t OL ) DACOS ( DL ) 

c ; ..... 

L1=DS0RM AlbAft**2 + 31BAR*<=2 ) 

2.L2 = DSQR T ( A2 8AR* *2 + B2.D AR**2 ) 

_____ DL = DSQR T I ( B2BAR-B 1 b AR ) * * 2*;. (A 2_t5AR-_A 1BA *2.1 

BE T 1 =DAR COS ( I L 1 **2 »DL **?-L 2 **2 ) / t 2 . «L 1 *0 1 ) ) 

B E T 2 = D A RC 0 S ( ( L 2 * * 2 <• l >L * . * 2 - L I* *2 )/ ( 2 ... * L£ *D L. ) 1 

EPS 1= D ARCOS ( <LHl**2<-DL**2-LH2**2) /( 2.»L H 1»0 L) ) 

EPS2 = D ARCUS I <LH2**2+DL**2-LH1**2 )/J 2.*LH2*0L Li 

ETA1=J.1615927-BET1-EPS1 

______ ETA2=3.1615927-BFT2-EPS2 

L0_= 010 R n Ll»*2+L H I »»2-2 . *L HI »L 1 »DCDS ( BET fVEPS 1 ) ) 

L A .M l = DATAN2(B1BAR,A1BAR_J 

F I LAMl .LT .0.0 I L AM 1= 6. 2831 853 1»LAM 1 

L AM2= D AT AN2 ( B2 BAR , A2 BAR) 

IFILAM2.LT. 0.0) LAM 2= 6. 283 135~ 31» L AM2 

L AMO=LAM2‘-DARCUSIIL2**2+LO *^2 -LH 2 «*2 ) / C 2 . »L 2*L 0 ) ) 

S I G 1 = 3 . 1615927-BE Tl-EPS 1-L AMI-*- L AMO 

IFILAMO .GE. 3. 16159 270) SIG1=SI 3 1-6. 28318 5 31 

. S I G 2 = 3 . 1 6 1 5 9 2 7 - 8 6 T 2 - E P S 2 r L A M 0 i AM 2 

I F( LAMO.GE. 3. 1 615 92 70) L AM 0= LA MO- 6 . 2 83 1 8 53 1 

_NU = L AMO 

__lam=nu 

A OBAR = LO^OCUS ( NU ) 
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BOBAR= LO *DS I N ( NU ) 

PHI 1= S I G1 

PH 1 2= S I G2 

( F I C 0 NF . E 0. 2 . )_ „K£T URN 

CALL S'JSPEN 

_W8! TCj.LDiXtJiii-C.UHI I 

MR I TE I I OUT , 54 ) 

' MR I T E I I HUT » 55) _ 

SSL S2= OSIN (S IG1+S IG2 ) 

C 2 S 1 2= DCOS ( SI 12) / SSI S2 

C l S 1 2= OCHS IS IG1I/SSIS2 

S2 S 12= 0SINISIG2) /SS1 S2 

S l S i 2 = OSINI SIG iF/SSI S2 

KBT = l./I (C2S12*C2S12/KH1) MClSl’2*ClS12/KH2>) 

KBL = t./( (S2S12*S2S12/KHl ) ♦ I Vl S 12*SJ S L2Z.KH2.lj 

KSPKHl=(2.*KKS( 1 ) *KMl ) / I 2 • *KKS ( 1 ) ♦KHl ) 

KSPKH2- (2 . *KKS ( 1 ) <=<112 )/ 1 2 . *KKS ( 1 ) *K H2 ) ; 

K = I 2 . * KK_S ( li * K. 8 L ) / 1 2 . * KK S 1 1 ) ♦ K B L ) 

BCT "iOP=fiEf I *57.2953 

8 E T2 DE = BE T2 *5 7 .29 58 ~ 

EPS1DE=EPSI*57.295 8™ " " 

EPS20E = E°S2 *57.29 58 

_ETA10E-STA1*57.2958 

ET A_2D E^E T A 2 * 3.7 .2 > 58 

S If, 10fi = S 101*57. 295d 

_SIG20E = SIG2*57.29 5 8 

_LAM10E=LAM 1*57. 2958 

_L AM2DE = LAM2 *5 7 .29 58 

L AM OD E = L A MO * 5 7 . 2 9 58 

NU O EG= N U 15 7 .29 5 8 

_ LAM 0 P n“= L A M OP * 5 7 . 2 9 5 8 

_NUP0EG.= NUP*57. 2958 

_EPSPID=EP$PI*57.2958 

_ EPS P20=EPSP2*57. 2958 

_ MR I TE ( lOUT, 56).. LH 1 , L 9 2 , A1 3 AR , A 2 BAR , B1 BAR , 82 BAR, LI , L 2 . DL , LO . BET 1 PE. 

,2L AM 00 E , NUDE G , A OB A R , H 0 B A R , K ., K S ?_KH_l_ t K S PK H2 . T I , 0T I . THE DU , T lj E DJ. . 

3TUR.DTVC.LS1.LS2.0LP. L0P.L AMQPQ.NUPD E G , AP BA R t BP BAR .EPS P I f) ,£ PS.P2Q 

4,0T ,KH1 ,KH2 

_ABAR = AJBAR 

_ RBAR= BO BAR 

XP=X»A6 AR»DCOSI THE)-BBAR*DS INI THEI-1 APBAR+L TR+.O 1 )*Dr.1S 1 THEP ) + BP8A 

_1R *0 S I N I THEP ) 

_ l P = Z * A B AR *0 SJ N I T HE )_* B BAR*DC3S( THE )- ( APB4R *LT R *.0 1 ) *D S I N I THE P ) - BPBA 

.1R*0C0S( THEP) _ 

_CONF= 2. 

_ RETURN 

END 

__S UB P. OU T INE SUSPEN 

_ I MPL I CJ T_ REAL *8 I A -H.3-2 ) 

DO UBLE PRECISIO N l AM I OE ,L AM2 DE , L AMO 0E , NJ DEG . L AMO P 0, K S . 

JK S P K H 2 , L T k , M U ,M Uf) ,L HI j L H2 . L 1 , L 2 . LO , L A MI , LAM2 . L AH 0 ,N. g 

2 , MUDEG , MU DDE G.NAP . NNP , LAM, LAMOEG. LS 1 . L S2 . LOP. L AMOP , NUP. NUPOEG 

3, M,MP,IY,IYP,LT,LTD,LT0,NA ,NN,K,KSPKH1 ’ 
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4, MMA116), MA,J IYPI 16 1 , KBL, KBT , KH 1 , KH2 1 K PH L» <KSJ_B ) 

COMMON J,3T ,X,Z ,XP ,ZP, THE , THEPfXD, ZD , XPO.,_ZP D, T MED, T HERD, CAM, GAMP, 

1 A L P , AL P P , AM, AMP, OYPR, OYPAP,RH3, S , SP , D.,.OP MP , t Y , I YP,LT,LTQ,LTD, 

20CG ,C ,K,CA,CN,CM,CMQ, CC A(b , 16) , CCN18, 16),CC M (8,16) I CCMQ(8,16),CAP, 

3CNP,CMP,CM3P,CCAP18,16) ,CC N° ( 8 , 1 6 ) , CC MP{ 0 , 1 6 ) ,CCMQP<3,16) ,V,VP,GR, 

4R i.AAi&jJl , Op.f 3 ,.31 , f U 31 , EE-l-il. DXaIL, Q XP.QZP .J THE . O T HE P . A P» A R . XB AR . 

5ZRAR,AAM(8) ,AAMP(R),AALPF(16),AALPP£( If,) , I .IN,. I OUT » 

60 AOTHE , DBDT HE , DAOTHP, DDDTH° , A3AR jBRAF. , A , B ,CH I , CH l 0 , MU , MUO , P H I ,SPl , 

7PHIP.PHI1 , PH 1 2, BET 1, BET 2, EPS1.EPS 2, KSPX HI, < SPKH2, L TK, DL , 

9S IGl , SIC 2 ,c TA1 .FTA2.TENS, J AMP , STHE , C T HE, STllEP,CTHEP ,LHL ,LH2 , ... 

9A1BAR , A2BA1 , B l .JAR , B2? AR , BO BAR , AOBAR ,1. 0, H ,L 2 , L AMO , L AMI , LAM2 ,NU ,G 

c OMMj.)N_rg K , T H E QUj.TH t P.L ..OUii i_LAM, II, LlN.FjUTI l£D6 Q 

l, Tin 16) , SSP I ( 16) ,L SI ,L S2 ,L0° , LAM'JP,NUP,BPBAH ,FPSPl .EPS °2 .OLP.KS, .... 

2,MMA,MA,TTG(8),VVC(8) ,V0, I IYP , <BL ,X8T ,K-1 1, KM2, TTENSI 8) ,KKS, 

3C0NT120) ,AALP(16) ,AALPP(16) _ 

C . 

C __ STATEMENT FUNCTION 

„darc;$.(.olp 1. (Di.Pl 

c . . 

E P SP1 =D ARCQS ( (01 P**2*LS 1 *-*2-LS 2**2 )/ (2.*0LP*L SI) ) _ 

EPSP2=OARCOS( (DLP**2* LS2** 2-LS 1**21/ (2.*OLP*LS2) ) . 

iop=n5ijRT nnip/2. )**2+lsi**2-oip*lsi*d;ds(sp$pu ) 

LAMOP=DARCOS< ( (OLP/2. ) ♦ *2 *L0 P* * 2-LS 2 **2 ) / ( DL P *L OP j _L„ 

N.y P_=L A.UjIP-1. 5 J . UJ 3 .bl 

APBAR=LOP*OCOS (NJP ) 

BP6AR=L0P*0SI N ( NU P ) . 

IF (CUNF.FQ.O.O) RETURN .... 

XP = X»ADAR*CTHE-BilAR*ST HE-AP 8 AR *CT HF PtBP3_AR*STHF P t ( L TOf TEN S/K)_ 

l*0S I N ( CHI ) . . _ 

Z P = Z «■ A B AR *S THE *B0 AR »CTHE- APB AR »STHE P-BPB AR»C THE P-H LTO» TENS/K) 

1*0C0S(C.MI) _ .. __ 


RETURN 

END 

SUBROUTINE CROUT ( A, C, N ,... LD, 
LINEAR ALGEBRAIC EQUATIONS - C 

.01 MENsioN A(L'),1 ) , C 111 

ZERO = O.OOQ 

ZERO = 0. 

I E R SW = 0 __ 

I F ( DABS ( A ( l , l ) ) - E P S ) 9 0 , 5 , 5 


LD. E PS. 
- CROUT 


.LERSW j 


LS220 _ 

hzzzzizi_z.z:zls2 2 o_’ 

. LS2 20 


5 

IF(N-L) 90,10,15 


LS22J 10 

10 

cm = C< 1 >/A! 1,1 ) 


LS220 11 

15 

RFTURN 


LS220 12 

DO 20 1=2, N 


LS220 13 

20 

All, I ) = A( 1 , !)/A(l, 1) 


LS220 14 


DO 65 1=2. N 


LS220 15 


DO 65 J = 2» N 


LS220 16 


SUM = ZERO 


LS220 1 7 


IF! J - 1 ) 30,30,25 


LS220 13 

_25 

J I M= I - 1 


LS223 19 


GO TO 35 


LS220 20 

30 

JIN =J-1 


LS220 21 


35 DO 40 K. = L , J I N LS220 22 



40 

SUM =SUM*A< I , K ) *A ( K,J) 

LS22J 

23 



IF (J-I )45, 45,55 

LS22J 

24 


45 

At I ,J) = All , J) -SUM 

L S220 

25 



IFl J-I )65, 50,90 

LS220 

26 ... ... . ... 

.50 I El DABS (A (1 . I) l-EPS )9Q. 65. 65 , 

C 

50 

I F ( ABSt M I , I ) ) -EPS) 90 .65.55 

L S220 

23 


55 

I F I DABS (All, I ) ) - EPS) 90 ,60 , 60 



C 

55 

I F ( ABS (All,!)) - FPS >90,50,63 

LS220 

3 J 


60 

A I I , J ) = ( A ( 1 , J ) - SUM ) /A ( I ,1) 

LS220 

31 


65 

CONT INUE 

LS220 

32 



cm = cm / Aii ,i) 

L S220 

3 3 



DO 75 1=2, N 

LS220 

3 4 



SUM = AERO 

L S2 2 0 

3 5 



JIN =1-1 

L S2 20 

36 



DO 70 K = 1 , J I N 

LS220 

3 r 


70 

SUM = SUM«-A< I ,K)*CIK ) 

LS220 

3 3 



C(I ) = (C(1 ) - SUM) /All , I ) 

LS2 20 

3 9 


75 

CUNT INOE 

LS2 2 J 

4 0 



JIN = N— l 

L S220 

4 1 



DO 85 M = 1 , J I N 

LS220 

42 



SUM = ZERO 

LS220 

43 



L = J I N-M+l 

L S220 

44 



LL=LH 

LS220 

45 



DO 80 K=LL, N 

LS220 

'♦ 6 


80 

SUM = SUM ♦A(L.K)*CIK) 

... ..LS2 20 

4 7 


85 

C(L ) = C I L ) -SUM 

..LS22D 

4-3 



RETURN 

LS220 

49 


90 

I E K SW= l 

LS220 

50 



RETURN 

LS220 

5 1 



END 

LS220 

52 

SUBROUTINE PL TR AJ ( X , Y l , Y2 , Y3 , Y 4 , N, IX , IY1 , IT 2 , I Y 3, I Y4, HDR ) 

C 


THIS SUBROUTINE PLOTS UP TO FOUR CURVES ON THE SAME X SCALE. 

EACH 


C 


HAVING ITS OWN Y SCALE 



c 

C 


1 ■ lilil JJTi JJ K JLI ® jl^TTYrnS^flEVflPTFfiZl^RlfSnnni 

[ <; rOMPITr 


C 


Yl.Y2.Y3.Y4 THE Y (VERICAL) COORDINATES AT EACH X 



C 


N THE NUM3ER OF POINTS IN EACH OF THE ARRAYS 

ABOVE 


c 


IX THE CONTROL FOR X LA3ELLINS (N0T=0) 



C 


IYl.IY2.IY3, I Y4 THE CONTROL FOR Y LABELLING (0=UMIT THIS £ 

FOLLOWING 

I 

c 


HDR THE HEADER IDENTIFICATION FOP THE CURVFSI 



c 

DIMENSION XI 1) ,Y1 11 > f Y2 (1 1 .Y3IL ) .Y4I1 ),HDR( l ) 

DIMENSION VMXI4) , VM1( 4) ,VM2t 4) . VM31 4) . VM4I4) .VMS I 4, 5) , I D ( 2 7 ) 


E qu I V AL ENCE .. < VM S ( 1,1 )., V Ml ( 1 )). ( VMS ( 1 . 2 ) . VM 2 1 1 V_M S< 1 , 1 ) . VH3 ( 111 . 


O D 
W ju 


• VELQC1 • ,'TY M /« SEC !.,? 

' ACCRUE '. 'RAT I3N ', » M/SE'.'C»*2 

IDYNAMI • ,'C P RES 1 , 1 SURE ■,' N /M«»2 

• TENS l 3 * , 'N N ' 

*FL IGHT v, •_ PATH_i, 'ANGLE » ,« DEG 

'PITCH • , •ATTITU* , • OE DE'.'G 


1__ __ (VMS! 1,4) .VMM 1)L, (VHS( 1 ,5) . VMX(l) ) 

DIMENSION LABELS I 4, 12) 

DATA LABELS' / • ALT I Tl) • , • UE M* • . • ID **3 




I F ( IY1.LE.0) GO TO 200 


IF! YU I J.GT.VMK2) ) 

VM1 (2) 

= 

vim 




IFIYU I I.LT . VM l ( I) ) 

VMKI) 

S 

Yl( I ) 




IF ( IY2. LE.01 GO Til 

200 






1F( Y2( I 1 . G T . VM2 I 2 ) ) 

VM2I2) 

St 

Y2 I l I 




J Ft Y2 U.J...L T «V M2 ( 1 ) ) 

VM211) 

S 

Y 2 ( I ) 




IF ( IY3.L6.0) GO TO 

200 





, 

IFIY3I I ) . GT • VM3 ( 2 I ) 

V M 3 1 2 ) 

=: 

Y 3( I > 




I F ( Y3 ( I I.LT .VM3 ( 1 ) ) 

VM3I1I 

= 

Y 3 1 I I 




IF ( I Y4.LE .0) GO TO 

200 






IF(Y4( l ) .GT .V M4 ( 2 1 ) 

VM4 ( 2 ) 

s 

Y4( I ) 




I F t Y4( I I. LT.VM4I1 I ) 

V M4 1 1 I 

= 

ym n 





200 CONTINUE 

c. 

C DO SCALING AMD AXES 

IF( IX - 1)400, 310,320 

310 VMX(l) = VMX(1)*FP 1" 

VMX ( 2 ) = VMX(2)*FP 


320 CALL SCALEIVMX.XL, 2, l J „ 

I" CALL AXIStO. ,0. .LABELS ( 1 f~IX ) ,-24,XL.0. , VMX I 3 I , VMX I 4) ) 

I F l IX - 1 1 340, 330, 340 ; 1_ “ 

330 CONTINUE.. 

VMX I 3 1 _= VMX 1 3 1 * F G ‘ 

VMX ( 4 1 = VMXI 41 *FG 
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CONTINUE 

X ! N ♦ l J_= V M X ( 3 ) 
XIN+2) = VMXI4) 


I F C I V 1 - 1)400,351,352 

VMI 1 11 = VMH1J»FP 

VMI (2) = VM1(2)*FP 

CALL SCALE! VMI, VL, 2, l L 

CALL AXISU. ,0. .LAPELS (1 . IV1).»24.YL. 

if! I Y l - 1) 354, 353,354 

CONTINUE _ 

_VM 11 31 = V M 1 1 3 I * F G 

VMI (A) = VM 1 l 4 ) *FG 

CONTINUE _ 

Y1INU) = VMI (3) 

Y 1 I N«-2 ) = VM1<4) 


JilJ.3 > J_ 


_I F ( I Y 2 - I I 400. 361. 362 

V M 2 ( l ) =_VM2(l)4-FP 

VM2 ( 2 ) = VM2(2)*FP 

CALL SCALE! VM2, YL , 2, IJ 

CALL AXIS! OA2.0., LAB EJL SCI. I Y 2 ) , *24, YL , 90 . . VM2 ! 3 ) , VM 21 4) ) 

IF! IY2 - l) 364.363.36 4 

CON TINUE 

VM2I 3) = VM2 ( 3 I +F G 

VM2(4) = VM2! 4) *F G __ 

.CONTINUE __ 

Y2 ( NH ) ^ VM2 ( 3 ) 

_Y 2(N»2). = VM2 C 4) 

~TfTTy~3 - 1) 400,371.3 72 

V M 3 ( l ) VM 3! I)*FP 

VM3I2) = VM3 ( 2 ].*FP 

CALL SCALE! VM3.YL.2.I ) 

CALL AX IS ( 3 A3. J.. LABE LS 1 I. IY3).»24. YL.90. ,VM3! 3) .VM314 ) ) 

IF I I Y3 - 1)374,373,374 

CONTINUE _ .... ! 

V M 3 ( 3 ) _=_YM3(3)*FG 

VM314) =_VM3!4)*FG 

CONTINUE 

_ Y 3 { N+l ) ..=. V M3 ( 3 ) 


IFHY4 - i jAoo.sai.jin. 

38 1 VM4I 1 > =_.VM4( 1 ) *FP 

_VM4 ( 2 ) = VM4(2)*_FP 

382 CALL SCALE! VI- 4, Y L, 2,1 ) 

CALL AXI SI 3A4.6T, LABEL! 

IF I I Y4 - 1 )3 84. 383 . 384 

38 3 C 0 NT I N U F. 

VM4 ( 3 ) = VM 4! 31* FG 

VM4 (4 ) = VM 4 (4 ) *FG 

_384. CONTINUE. 

Y4INUI = V M 4 ( 3 ) 
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CALL L INEI X, Y4,N, l, J, LTYPE l 41 ) 
CUMTI MUE 


OP AW LABELLING _ 

VX = XL - .5 

T = V X -2.7 

VY = YL - HT-HT _ 

IF (IY1 >550,550,510 

CONTINUE _ _ 

CALL SYMBOL (T,VY, HT ,L AO EL Sj l , IYiLi0.j2 
CALL SYMBOL (VX, VY+HT2 , HT , L TY PE UJ , 3. , - 

V Y = VY ~ UT-HT 

1>(IY2 )550", 550,520. _ _ 

.CONTINUE 

CALL SYMBOL I T , V Y, HT , L AU EL S ( 1,1 Y2),0. ,2 
.CALL SYMBOL (VX,VY+!IT2 ,HT, L TYPE.! 2),. 0_.,~ 

VY = VY - HT-MT _ 

JFJ IY3) 550, 550,530 

CONTINUE .... ... .. 

CALL S YMT OL I T , V Y , H T , L Afi EL S ( l , I Y3i j,0. l 2 
CALL SYMBOL l VX, VY *-HT2 , HT,.L TYPE ( 3)., 0,.,.- 

VY = VY - HT-HT _ 

IF< IY4) 550, 550, 540 

CO N T I NUE 

CAL L ~S V MBOL ( T , VY , HT , L ABEL S 1 1 , I Y4J.,0 . , 2 
CALL SYMBOL I VX, VY+HT2 ,MT,L TY PE ( 4> .0. 
CONTINUE 


CALL PLOT I XL+4,,-1. 5,-3) _ 

.CON 1 1 NU E 

.RETURN 

END 

SUBROUT INF OINSMUl , P , RHO . CS I 

DOUBLE PRECISION Zl, P, RHO, CS 

_0 1 MENS ION „ HB 1 22 ) , _P B l 22 » , _ JB ( 2 2 1 , AI22> . 
DATA HB/Q., 36 08 9.2 39, 65 616. 798, 1049 86, 83 i 


B ( 22) 

154199 « 48 , 170603 . 67 , 
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1 200 131. 23,259186. 35. 29 H 5 3. 4 0 , 3 23 002 ♦ 7 5, 3 54 75 3. 59,386406.39 , 

' 2 4 80 731.04,8 120 4 6. { 6 ,' 54 3 21 5.40,805263 .46, 7 2 9 24 3 .9 1, 9 39894.7 5. 

3 1234619. 4 . 1520799.4,1793726.4,2068 776 .6/ 

1_ DATA P.3/ 21 16.217, 472. 67922,.i 14.34505, 18. 128852,2. 3162994, 

1_ . 1. 2 322512 ,. 380321 73 ,.021672818 , .0034331482^ .00062912953t 

2 . .00 31 5 3 55,19 36, .(30 O.QS 2666807. .000 0 10 57 15 8? .. 0000 0771 5 70 71 ■ 

3 .0000 05 3324672 , .0000 0 3 6196 1 39 , .000001.45.3.71, .39343987E-6, 

__ 4 . 941 7o667E- 7, .22884 174E-7 ,. 720689369-8, .248912648-9/ . 

DATA TQ/298. 15, 2*216. 65, 228. 65, 2*270. 65, 252. 65, 2*180. 65, 210. 65, 

1 260.65 , 360. 65 ,560 .65, 1 U 0.65 , 1210.65, 1350,65, 1550.65, 18 30.65. 

2. 2 16 0. 65, 2420. 66,2 590. 65,2700.65/ _ 

OAT A A/- . 6 3_7 6 5 356 r-6 1.0 . L , J 4 6 3 775E- 5, .37 J 251 60F 

1 -. 2252 3 5541: - 5 , - .4 8266491 E - 5 ,0 . , .5 2 1 4 1 4 J8H-5 , . 74757236E-5, 

2 . . 1212076904, . I 7623 20 1 E-4 , . 4994 1 99 7F - 5 , . 2839653 7E-5 , _ 

_ 3 . 18835746E-5, .12041 172C-5, .85314774E-6, .6H63ft7E-6, _ 

4 . 4 2043419E-6,.25268S39E-6»2*.l5723l5E-6/ . . . ._ 

DATA 0/ 5.25 58 86,- .48063 1023-4, - 34. 163232 12.201 l 70 ,-.334 73567F.-4 , 

1 17 ,_0 015 2 7_, a .54 0804. -.8 76 41 136E -4. - 11.05 5 22 6, - 6. 6ji 2790 1 . 

_ 2 -3. ’2961763, -1.6390358, -2. 1704464, -3. 2456970,-4.6156949, 

3 -6. 4033863,-/. 3733154, -9. 3039039^17. 4627,- 17,025193,, 

4 _ _ 2*-25. 562133/ 

~Z" Z ” Z 1*3.28084 I ’ TZ 

H * 20 3755 3 1 . F 0*Z/( 2085553 1.E0 till I 

_DQ 1 I = 2, 22 

I Ft H-H8 ( I )> 2 , 1 , 1 

1 CONTINUE _ 


_ I = 23 _ _ 

2_ . j_ _ = i -i 

DH = H - HB ( j i 

TEMP = 1. 4 AU)*DH 

T = TO ( I)*TEMP 

IF(A(.m3,6,3 

6 TEMP = EXP(B(1)*DHJ 

GO TO 4 

3 TEMP - T E MP**B ( 1 ) Z 

4 P = P5 U ) *T EMP»47. 88025 

__ RHO = .003483647* P/T 

... CS = 20 . 046 79* S03 T ( T» 

_ RETURN 

END 


Page 76 
GER 15853 


- iO\ *N-A 1U/ 

: ENGINEERING PROCEDURE S-017 



lo Lagrangian Dynamics , D „ A, Wells, Schaum Publishing 
Company, New York, No Y„ , 1967 

2 o Classical Mechanics, H. Goldstein, Addison-Wesley- 
Cambridge, Mass., 1950 

3) Introduction to Numerical Analysis, F.B. Hildegrand, 
McGraw-Hill Inc, 1958 

4) Mathematical Handbook for Scientists and Engineers, 

G„ A , Korn, McGraw-Hill Inc, 1961 

5) Mathematical Model and Computer Program for a 6 D.O.F. 
Body Connected to a 5 D.O.F. Body by an Elastic 
Tether, George R ( Doyle Jr. GER-14958, August 1970 




